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Convergencé of finite element solutions in 


viscous flow problems (") “ 


Convergência do método dos elementos finitos 


aplicado aos escoamentos viscosos 


RESUMO 


Resumem-sa primeiramenta os resultados de uma 
teoria dos métodos variacionais cuja utilização permite 
a análise da convergência do método dos elementos 
finitos aplicado ao escoamento de um fluido viscoso 
regido pelas equações de Navier-Stokes. O método em- 
pregado requer a decomposição do fluído em dois corpos 
mais simples: um fluído viscoso sem inércia e um 
fluído com inércia mas sem viscosidade. Apresenta-se 
para este último um novo princípio de minimo. 


INTRODUCTION 


The problem of the convergence of finite element 
approximations has been considared by the author in 
several papers, e.g. Reference 1, connected with struc- 
tural analysis. 

in a more recent paper,? a theory of variational 
methods is presented, general enough to cover any non 
“linear, non-structural cases for which a minimum prin- 
ciple is known to exist. 


E. R. DE ARANTES E OLIVEIRA (2) 
Departament of Civil Engineering, 
Technical University of Lisbon, 


Lisbon, Portugal 


ABSTRACT 


Results of a theory of variational problems are first 
summarized. The theory is then applied to the analisys 
of convergence of the finite element method applied 
to the flow of a viscous fluid governed by Navier-Stokes 
equations. The approach requires the decomposition of 
tho body into two simpler ones: a viscous non-inertial 
fluid and an inertial non-viscous one, the mechanical 
behaviour of which is governed by minimum principles, 
A new minimum principle is presented for the inertial 
non-viscous fluid. 


lt will be shown along the present paper how 
such theory can be applied to the discussion of the 
convergenca of finite element approximations in viscous 
flow problems governed by Navier-Stokes equation. 

The trouble with Navier-Stokes equation is that 
it cannot be associated with a variational principle in- 
volving only the velocities, although a variational prin- 
ciplo can be found for the system formed by Navier- 
-Stokes equation and its adjoint equation.? Such a prin- 
ciple cannot be of much help for proving convergence, 


[1 The research reported in this paper was supported by Instituto de Alta Cultura, Lisboa (Projecto TLE/4). 


(2) Professor. 


(*) (Publicado pela primeira vez no «international Journal for Numerical Methods in Enginesring», Vol, 3, 739 


— 163, 1975). 


[*) First Published in tha alnternational Journal for Numerical Methods in Engineeringa, Vol. 9, 739 — 763, 1975). 
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however, because it is not a minimum principle as our 
theory requires. 

Fortunately, however, the theory can still be ap- 
plied to complex models not governed by global mini- 
mum principles if such models can be shown to result 
from the combination of simpler models each of which 
is indeed governed by a minimum principle. 

This is precisely what happens with a fluid flo- 
wing in accordance with Navier-Stokes equation, as 
we can consider it as resulting from the association 
in parallel of two simpler bodies: one with viscosity 
and no mass (viscous non-inertial flow), tha other with 
mass but no viscosity (inertial non-viscous flow). Such 
association is schematically represented in Figure 1. 

The minimum principles of the theory of viscosity 
permits discussion of convergence in the first case. 
A principle closely connected to Hamilton's principle 
of Mechanics permits discussion of convergence in the 
second case. 

The main results of the theory will be summarized 
in the next section. lts results will be applied to the 
viscous non-inertial flow case in the section “Viscous 
non-inertial flow and to the inertial non-viscous flow 
case in the section 'Inertial non-viscous flow', A theo- 
retical investigation of how the results obtained in the 
sections Viscous non-inertial flow' and 'Inertial non- 
-viscous flow can be used for the discussion of con- 
vergenca in the general case of the inertial viscous flow 
will finally be presented in the section Viscous iner- 
tial flow”. 


RESULTS OF THE THEORY OF VARIATIONAL 
METHODS 


Let X be a Banach space and let q be a continuo- 
us functional on X which is assumed to be minimized 
on a given subset( c Xby an element s EC. 

Our basic variational problem consists in the de- 
termination of such a minimizer, s, which is called 
the exact solution, 

s is assumed to be a critical point of y on €. The 
first and second Fréchet derivatives of qy are supposed 
lo exist at s and the second derivativa of q to be a 
positive definita operator. 


Non -inertial 


Figure | 
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A projection operator B, which may be called an 
inter polation operator, is considered which makes any 
element xe X correspond to a certain element x” of 
a subspace À Cd. 

Let C' o X' be the set of all the B-images of the 
elements of C. It is assumed that an element se € 
exists which minimizes pq on C'.s, is called the appro- 
ximata solution. 

The approximation theorem”? states that a metric 
can be introduced in X such that, if infinitesimals of 
higher order are neglected, 


ds, 54) < dis, Si) + (|ôpl +lô49|) (1) 


where s, is any element of C such that 


Ss, = B(s,) (2) 

and 
dp = «p(s)—qp(s) (3) 
duP == P(sa)— p(s,) (4) 


C' does not generally belong to C. If it does, as 
in tho case of Ritz's method, operator B is called con- 
forming. 

In the conforming case, s, can be chosen to 
coincide with s, and inequality (1) becomes simply 


d?(s, 5.) < |ô,9) (5) 


Now, it can be shcwn? that, as se É. 
dp = dº(s,s). (6) 


if infinitesimals of higher order are neglected, so that 


t+ yields d(s.s.) < dis.) (7) 


Inequality (7) means that, in Ritz's method, the 
approximation error is bounded by the interpolation 
error, 

In the more general case of the finite element 
method, C' is not a subset of C, so that the use of 
inequality (1) is required in the discussion of conver- 
gence. 


VISCOUS NON-INERTIAL FLOW 


The equations of the linear theory of viscosity may 
be obtained from those of the linear theory of elasti- 
city just by replacing strains and displacements by rates 
of strain and velocities. We have then 


Ci P;tpf,=0 (8) 
et 300, = Hb, ;+D;i) (9) 
cj; = 2uê;, (10) 

p= —Kké (11) 
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The stress and the strain rate tensors are decom- 
posed into their isotropic and deviatoric terms accor- 
ding to 


0; = Cij—P Oi (12) 


ê, = êp+IdO,, (13) 


The domain A is assumed subdivided into sub- 
domains 4. 

The boundary of each subdomain Aº* is denoted 
by Bº and the external boundary of 4 by B (subdivided 
into B, and B,). The internal boundary of à is denoted 
by |. 

The boundary conditions on B are 


3; =Gn—pn,=6, on B, (14) 
whera o Í denotes the component j of the traction vector. 


On the internal boundary we assume 
o; +o, =0 (16) 


i (17) 


where indices + and — denote each side of the inter- 
nal boundary. The first equation means that no external 
forces are distributed on [ and the second that the 
velocities are continuous across the subdomain boun- 
daries. 

A pair of stress and strain rate fields connected 
by tha stress-strain rate equations (10) and (11) will 
be called a viscous field (or simply a field along the 
present section). 

A velocity field connected with the strain-rate field 
through the strain rate- velocity equations can of course 
be associated to each viscous field. Such a velocity 
field, however, presents 6 degrees-of-freedom wich can 
be eliminated only if adequate velocity boundary con- 
ditions are given on B,. 

The elements of space X are in the present case 
all the viscous fields on 4. The norm of an arbitrary 
element xe X is given by 


xi = [=] istydat) = Je. (és + é dar] 
(18) 


Two minimum theorems appear in the theory of 
viscosity which correspond to the minimum total poten- 
tial and total complementary theorems of the theory 
of elasticity.! 

The complementary energy approach will not be 
considered in tha present paper. The functional q is 
therefore the functional involved in the first theorem, 
i.e, tho total energy dissipation rate 
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A de z Ef, tor; to Cj H — pé) dA“ | pf, dA! — 
R vº Ae 


O: am | (19) 
“Roo B, 


The variational principle is an Eulerian one, l.e., 
tho exac: solution makes the first Eulerian variation of 
py equal to zero on the set of all the viscous fields 
which satisfy equations (9), (15) and (17). 

The subsei C  - Xis thus the set of all such fields. 

Let us consider now a decomposition into finite 
elements cor.esponding to the subdivision into subdo- 
mains and let a set of points on the internal boundary 
bo selected as nodes. 

Let the elements of X' be fields which satisfy 
equations 

UV; = WU (20) 


1 a 
Ci > XijanDaN (21) 


within subdomain 4º. 

index N denotes a general node of the system 
and v;y the nodal values of functions V',. 

The elements of C' are the elements of X' which 
satisfy the bounda.y condition 


Din = CAN) (22) 


for nodes located on B,. 
The discrete equation which represents the condi- 
tions for qy being stationary on Cº Is 
a 
Mimnlim = Jin (04) 
where 


MiMnN = 3 
u A” 


[LUX im esiN — (4 UF RX rim Ls jn] |dAº 


(24) 


fin => f UNPÍ, dar+ [o , Wo, apr (25) 


u 


Considering the results obtained for the elastic 
case in a previous paper,? it is easy to conclude that, 
if x is an arbitrary element of X, x is the B-projection 
of x in X' and | is the largest diameter in the whole 
set of subdomains A*;.the variation of q from x to xl ji.e., 


dp = plx)— (x) (26) 

is of the order of |"* ! if 

(i) n is the degree of the highest degree complete 
polynomial with arbitrary coefficients contained in any 
of tho expressions of the right-hand sides of (20) end 
(21). 

(ii) all the derivatives of the velocities and strain 
rates associated with x, with order n + 1 or less, are 
bounded within each subdomain Aº. 
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The distance «l(s.s,) between the exact and the 
approximate solution is of the order of |!"' | if con- 
dition (i) is satisfied and condition (ii) is also satisfied 
both by the exact and approximate solutions* (whate- 
ver the value of |). 


In case of conformity, (5) holds true and, there- 


fore, dls,5,) is of the order of |"*!!= whenever con- 
dition (i) is satisfied and condition (ii) is also satisfied 
by the exact solution s. 


INERTIAL NON-VISCOUS FLOW 
Equations 


In absence of viscosity, the stress tensor g;, is isotro- 
pic, i.e., 
Oi =— — pô; (21) 


where p is the pressure. 
The equilibrium equations become therefore 


pótp; = pf (28) 


where /, denotes the body force density. 

These 3 equations involve 5 unknowns: the three 
velocities, the pressure and ths density p. 

The two equations which are missing are the con- 
tinuity equation (density-velocity equation) and the 
equation of state (pressure-density equation). 

The continuity equation can be derived from the 
condition 


d 
“todAj=0 (29 
dt? 
which expresses the invariance of mass. 
Às 
d (d/dij(dA) - du | 
— 4 Eua ps | —— — Y dA. 
griSA) RA dA =" EP dA = divy 
(50) 
(29) vields 
dp Ê 
= dA+pdivvdA =0 (31) 
and, therefore, 
d ais 
AR divv = 0 (52) 


di 


which is the Eulerian form of the continuity equation. 
Another form of (32) is 


Cp 
= Gy = Q [54] 
Ut 

Let us consider now the equation of state. 

We start by assuming that no mechanical energy 
is dissipated into heat, and writing that the variation 
of the internal energy contained within an element of 


mass dm is equal to the work done by the external 
forces, i.e. 


ó(e dm) = —p ó(dA) (34) 


where e is the internal density. 
By virtus of the invariance of mass (continuity), 
(34) vyields 


dm dl ss 
de dm = “pé | = “pol dm (35) 
ê, É! 
and, therefore, tho equation of state 
de pila 
= 40d 
Pp d(1/p) | 
or” 
ds 36b 
Rm 
f Í dp (350D) 


Density e is admitted to be such a function of | 
such that the continuity of p implies the continuity of p. 

lf tha domain A is decomposed into subdomains 
Aº, the exact solution is supposed to satisfy the equati- 
ons (28), (32) and (36a) within each subdomain As, 
the external boundary conditions 


) E 
Ea [| on B, (57) 
P É 

Um — a om B, [38] 


(394) 
L, = Um (40) 
on 1 


As the continuity of p implies the continuity of p, 
equation (39a) implies 


Wla à 


which is analogous to (37). 

Un denotes the normal component of the velocity 
at a given point of the boundary. Equation (40) supposes 
that the positive sense of the unit normal vector at each 
point on [ is defined (from — to +, for instance). 


(39b) 


A set of a velocity field, a pressure field and a 
density field will be called an inertial field (or simply 
a field, in the present section). 

A displacement fisld with components uu, such 
that 


à dt = U; (41) 


* lt can be demonstrated that if condition (1) holds and the so-called “patch test” is passed by the system of elements, 


condition (1) will be satisfed by s.. 
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can of course be associated to each inertial field. 

An inertial field will be termed continuous if equa- 
tions (32) and (36) are satisfied within each subdomain 
àº and equations (38) and (40) are satisfisd respec- 
tively on B, and on T. 


Variational principle 

Different variational principles have been presen- 
ted by several authors*5 for the case of non-viscous 
flow. 

Such principles can be classified into Eulerian or 
Lagrangian, according to which kind of variation of the 
functional is made equal to zero, 

This classification of the principles into Lagrangian 
and Eulerian cannot depend of course on the functional 
baing expressed in a Lagrangian or Eularian form ba- 
cause although some expressions are more naturally 
written in one of these forms than in the other, any 
expression which is written in Lagrangian form can also 
bo written in Eulerian and vise versa. 

Now, the problem of integrating equations (28), 
(32) and (36) over a fixed domain A is an Eulerian 
problem in the sense that it considers at each instant 
the portion of the fluid actually contained in A and 
does not follow in time the motion of a definite mass. 

lt seems thus in accordance with the nature of 


such a problem that the functional to be made statio- 
nary has a natural Eulerian expression, involving inte- 
grals on 4 of Eulerian functions. 

This does not mean, however, that an Eulerian 
principle, like Lin's principle,” has to be used, because, 
as pointed out above, what makes a principle Eulerian or 
Lagrangian is the kind of variation which is made equal 
to zero and not the expression of the functional. In 
other words, a Lagrangian principle can be associated 
to a functional whose natural expression is Eulerian. 

A Lagrangian principle can will be used therefore 
which has the advantage over Lin's principle of not 
involving Lagrangian co-ordinates as constraints. 

The functional is 


o =| (r- V) di (42) 
where 
T = 1 | puiv; dA (43) 
E [a ea 
V=E-) Í fix;p dAº+ [ [3 xp dB (44) 
e das Bb 


T is the kinetic energy and WV performs the role 
of the total potential energy, so that the principle is 
closely connected to the classical Hamilton principle. 
E is the internal energy. x, denotes x,n, 

Our variational principle states that the exact solu- 
tion of equations (28), (32), (36a) and (41) with initial 
and final conditions 
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uAt,) — uAt,) (45) 
uáto) = ut») (46) 
plti) = pl) e 

(48) 


lts) — Mts) 


and boundary conditions (37) — (40) makes q stationary 
or, more precisely, makes tha Lagrangian variation of 
y equal to zero on the set of continuous fields on 4 
which satisfies (32), (36a), (38), (40), (41) and (45)-— 
— (48), 

The fact that variation is Lagrangian and continui- 
ty is preserved implies 


do. = Sá no (49) 
&p dA?) = k within Af (50) 
ôv, = O on B, (51) 
dv; =00 = 00, of F (52) 
Also, by virtue of (45) — (48), 
dudt,) = O (53) 
dudto) = O (54) 
dm) = 0 (55) 
dp(to) = O (56) 


The Lagrangian variation of the kinetic energy is 


AT = Iz | Apr; dA") = 5º |] 300, Ô(p dA) + 
e dar Nº 


É 


o e a! 


di 


a 


pe,ôóv dA = sy | pr; du; dA'= 
“ xe 


i A 


d 
- GE 


| .d 
pr; du, dA — 5 dus (pri dA”) (57) 
” A 


a o 


By virtue of (29), 


d d 
dr": dAº) = Vig? dA“) + vp dA* = tp dAS, (58) 


there results 


d 
DI = di ba | podu dA — 5 | óôdui,p dA: (59) 
e “AF 


E A” 


143 


On the other hand, using (36b), we obtain 


dE = 0 | cp dA = | dep dA +| có(p dA") - 
A“ qe 


“o 


E ce. nm 
| dep dAº = — dp pdA' = op dA: [6()) 
.p 


dy dia Cp 


and, as (32) vields 


op+pdivôu = 0, (61) 
thare results 
Ô | ep dA* = | pdivóu dA = — 
4” à hs 
= | div (p ôu) das + | grad p. du dAº 
“Ao Ne 
=—| pôu dAS+ | grad p. ôu dA (62) 
a 4a 


As, by virtue of (49, (51), (52), and also (45) 
and (46), 


du, = 0 on B, (63) 
du, =ô0u, =ôu, on T (64) 
wa can write 
T | póu,dB' = | (p-p*)ôu,dr + | pôu, dB 
cv “HR I “H, à 
(65) 
and, therefore, 
DE = O Z| ep dA = | p; ou, dA 
Fi 4º o ad 
+ | (p —p )Jôu,dT — | pôu, dB (00) 
“| H 
Considering that dx, = du; 
O | fixpdA* = | fi du,p dA' + 
uy a” 
+ | tu; ó(p dA) = | hôup dA (67) 
o Je 
* p a p 
O | [-|x,p dB = | du,p dB+ 
nb Bob 
a + l h 'p) 
+| u, lp dB) = | ôu,p dB (68) 
“ LE P, “ni, H 
Introducing (59), (66), (67) and (68) into 
Op — | (0T-ôoVIdr. ht] 


Fi 
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we obtain finally 


dial ES 
dp = | | ) | pr, Ôu, 48] = 
Na da 


pi 


= > | (pr, + D;— pf) Ou, da“ 
vo vv AE 


es 


Pp IPs 
+ — | OU) dB, = tp E VOU, di | di 
“nm, LÊ pt 


(70)) 
By virtua of (53) and (54), 
| 5 | pr;ôu; dA] di = 
NET es 
; é — | = 
= | Pr; OU; as = () (71) 
e Lad Mi 
and thus 
Mp = | [5 | (pf, — pê;— pPaldu, dA 4 
bi O 


Ji | 20)| du, dB- | (pp ud dr (72) 
JH, É! b! | [ 


Equations (28), (37) and (39a) are therefore the 
necessary and sufficient conditions for functional q 
becoming stationary. 

In what conditions will the variational principle 
which has just been introduced be a minimum principle? 

A similar discussion has been made for Hamil- 
ton's principle (see, for instance Reference 7, pp. 651- 
— 667) and could practically be repeated for the pre- 
sent case. 

The essential result of this discussion is that the 
principle is a minimum principle if t, — t, is sufficien- 
tly small. 


Indeed, we obtain from (72) 


| Ej — E du; + "a su [o dA + 
ti e Ar P 


+[ (8 ÓU, dB | (0p' —-ôp )óu, dT ai (73) 
H, r 


Il 


(520), 


p 


where s denotes the exact solution. 
But 


ad Es : 
— dt; ÔU; = — Frade du;) + dt; ÔU, (74) 
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and thus. by virtue of (53), (54) and (29). 


mf 


-- I2 
(0“9), = | ba ) pôóv, óv, dAº di+ | > = 
bj E Ar E E 


sm | (23 du;p dA" + 
' p 


o) e 


| to aB- | op" —ôp" You, dr] dt (75) 


+] 5!” 
H 


po NV 


Considering (63) and (64), we can write 


| (2 ôusp dB + | 5) -5(2 | ouso dr — 
“B, p IL Pp p 


a | o(2) ôu,p dBº (76) 
é Sar iy 


Bv virtue of (36b) and (61), 


(8) | 5 Er e] ; 
Poe] = polpsl=ps TEsa| p= 
E del. = sd*e õp 
2 + 08,3 ni [p+o | p 


= — [AA divôu = [p-0? div du (77) 
| 3” | dp 
so that 


| (8) du,p dBº = | [p- ng div ôu óu, dBº 
- AP o dp 
(78) 


On the other hand, 


p (23) = dp;tpadiv du = E dp 
p dp 


+p.div du 


d 
= fps? div du +p. div du (79) 
so th 


ns ) | dp .. | o fai 
- vi] óu,p dAº = —— div du| du, dA*— 
JA | p | Ae “do , . 


| pa div ôu du, dAº 
A” 


= | [02 iv ou su] das | 
se | dp F As 


— | padivóu du, dAº 


A” 


pa rdiy du)” dAº — 
dp 
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e OD, o 
ei —— div du du, aBe— | ——[(div du)” dAº — 
Br "do A “do 
—| p; div ôu ôu dA (80) 
a E fed 


Introducing into (75), there results : 


(OP), = | E | pôr, ôóv dA“ dr+ | re 
[| Fi A” 


“|, 


dom 


R Dad o | 
=> 1 | — grad podu=poo div su div du dA 
Nº b 


U 


+] 
Hr 


It will be shewn now that the first term is of lesser 
order than the remaining ones so that, as such a first 
term is always positive, (0"p), is also positive if t, 
— t, is sufficiently small. 

Indeed, 


dp. ã | 
p ôu,— pa su, div ou am] di (81) 
bt 


dult) = duty) +(t— to) duto) + à 


A | 
“DL súgr) (5% 


where t, < t' < tand thus, by virtue of (49) and (53), 


, Us o o 
On the other hand, 
OUÃt) = d0dty)+(t— t,)dUÃt) (Na) 


where t, & t" S 1. 
Multiplying (84) by (t-t,) and subtracting from 


dulo) = (t— ty) d0AD + — EIS OBA) CAD DO (NS) 


Assuming that the fields of the class on which 
the functional is made stationary have bounded acce- 
lerations, the last term in the righthand side of (85) 
is of order higher than the order of the first so that 
we can write 


OU; — (4 — 1, )O(dr,) (86) 


Assuming that the unit volume expansion is also 
bounded for such fields, 


ô(divu) = O(ó ju!) (87) 
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Ui PM 


and then 
o ts 
| 5 | Õr, OV, dA! > | > | — 
“tj E A ou 
( dp | . f” : é 
] po, tdiv ou)" dAC4 | tdi du) du, dBº | di 
co Gp Rm 
(Ns) 
lt becomes clear that 


(Op), > O (NU) 


whenever ||dv| > Oand (t, — t,) is sufficiently small. 

Assume now that such a small interval of time in- 
creases and reaches a critical value above which(9"9). 
ceases to be positive. | 

For this critical value (0). = 0, which means 
that several different motions are possible and that the 
motion becomes unstable, 

In other words, if the location of the fluid parti- 
cles is prescribed for t = t, and t = t,, the motion is sta- 
blo if t, —t, is sufficiently small. lf t, — t, increases, 
the motion may become unstable. The important point 
is that (dq), can become negative only in case of ins- 
tability. 

lt must be added that, although this conclusion 
seems to be quite reasonable, an investigation on the 
connection between the sign of (5), and instahility, 
in particular: cases, is certainly lacking and should be 
made in the near future. 


Discrete equations 


Let us consider the body decomposed into Eulerian 
finite elements corresponding to the subdomains As. 

Although our analysis could be adapted to elements 
with other shapes, it is assumed in the sequel that such 
elements are tetrahedral. 

Two kinds of nodes are considered: the vertices 
and points located on the edges of the tetrahed;a. The 
generalized displacements and velocities are ascribed to 
the ve:tices. The generalized densities are ascribed both 
to tho vertices and to the mid-edge nodes (Figure 2). 


Um ' Vs PNN 


Eh HP 


PLM 
Figure 2 
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un Un and Pwx denote the values of magnitudes 

u;. vt; and p at node N. Pp; w denotes the value of the 

density p at the mid-edge node located between the 
vertices L and M. 


We assume 


Axo) = Wade) (90u] 
poe) = piu uti) (91) 
plxoo) = ptx Oyster yr, ut) (92) 
within As, 
Function 1/,(x) is assumed such that 
WXy) = dm (93) 
and 
x, =— YX)X;n (94) 
As, on the other hand, (41) is supposed to hold, 
OU, = WNÔU;N (90b) 
Function q), ylXx) is such that 
PimiXas) = ÔLM.RS (95) 


whero d,m ns is assumed equal to one if L coincides 
with R and M with S, and to zero in any other case. 

Magnitudes 7,« are not nodal values of [' ||, be- 
cause the product 1//,14'y does not take values one and 
zero at the mid-edge nodes, 

On the other hand, magnitudes rr, are the only 
nodal magnitudes which are denoted with index e. This 
is duo to the fact that the nodal continuity of u,.r, 
and ;1, and not of p'js, is postulated. 

The generalized form of the continuity condition 
(29) is admitted to be 


pm 
F, 
as 


Ud 


ei | | 
| Wi qo dA") = 0 (96) 


Proceeding as in (27), we obtain from (96) 


5] UP RAS e a p dir ] dA = 0) (97) 
PAR di 
and 
varia TOP a no | 
l sr ti 


Introducing (902) and (91) into (98) we obtain 
tho discrete or generalized continuity equation 
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dimpoP ro + Dimpanfpqlin = O (99) 
where 
ULmpo = z) VibuPro dA* (100) 
e “Ar 
human = |. WiWidorobs), da: (101) 
E Ne 


We observe that the number of non-trivigl indepen- 
dent equations in (99) equals the number of nodal den- 
sities ppo » i.e. the total number of vertices and mid- 
-adge nodes in tha tetrahedron. 

Indeed, as the coefficients u, mpq and b;;mpon are 
symmetric in P and Q, equation (99) does not change 
if P and Q are interchanged. On the other hand, all the 
equations (99) corresponding to two nodes P and Q, 
such that no element exists to which P and Q simulta- 
neously belong, are trivial, because the product pq . 
and thus the coefficients a, mp and Db;mpqn are then 
equal to zero. ) 

As equations (91) and (92) are not compatible 
with the equation of state (36a), it follows that this 
equation is not satisfied within 4º. 


A generalized equation of state can be established 
if a Lagrangian variation 


dE" = —| po(dA') (100) 
A” 


of the internal energy of a general element e, corres- 
ponding to a Lagrangian variation of the density of the 
form 


dp = PimliX)OPim (103) 
where dp, m denotes a Lagrangian variation of p at node 
LM, is admitted to be the same, whether it is computed 
with the distribution of pressure ylelded by the equation 
of state (36b) or with the distribution defined by (92). 


In the first case, we have, by virtue of (36a), (50) 
and (102), 


de 


de : 
JE = — 2 — (dA) = dpp dA" = 
í É I. b dp' ( ) " dp Pp 

de á AME w 


In the second case, by virtus of (91), (96) and 
(102), 


0E* = -[ vibyp o dA riu = 
nº 


=[ Viym dp dA'Tiy = | Viu Pro dA ÔPpoTiM 
se A 
(105) 
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Comparison of (104) and (105) vields 


| Wibm Pro dA'nim ôPrq oe 
A” 


de . jolie | 
IR dp" Pro dA“ ôdppo (106) 
and, as ôdppoq is arbitrary, 
SimpoTim = hpolp) (107) 
where 
£LmpQ — | Viyupro dA (108) 
Ae 
pe de na 
as Gp 


Equation (107), which contains indeed as many 
equations as unknowns 7, - is the generalized equa- 
tion of state for element e. 

An important point is that 5p given by (103) is 
not always a Lagrangian variation of P on the set of 
functions defined by (91). 

Thc expression of such variation is indeed, by vir 
tue of (90b) 


| 


PRsiX) dE Past PRs AX) OU; Pas 
pRÁX) de Prs + PRsOOVNOOPRs ÔUin 


ôp 


| 


(110) 


where d,p;m denotes an Eulerian variation of P at the 
node RS. 

The Lagrangian variation of p at node LM is thus 
given by 


dpim = PRÁLM) depas + Prs LM ALM pas ÔU;N 
(111) 


which, by virtue of (95), can be transformed into 


dpim = dePim + Pre LM ALM Ipes dus (12) 


Introducing (112) into (103), there results 


dp = PiulX) depim + Piu )PRs ALM A LM pes ÔUt 
(113) 


Comparing (113) with (110), there results that, 
in order that the two expressions of 5/ be the same, it 
is necessary that 


Prs OU Nd) = Pipas ALMADA, (114) 


14% 


This means that each function Pes AX)yX) 
must have the same distribution on Aº as the density p 
(seo (91)). 

Such a condition can easily be satisfisd as we are 
going to show. 

We observo that, considering the number of verti- 
ces and mid-edge nodes, magnitude r, can be given 
an arbitrary linear distribution, while p and p/p can each 
be given an arbitrary quadratic distribution, within 4º. 
Functions prAxX) will therefore be quadratic, functions 

W (x) will be linear and functions v; «dxvios will be 
quadratic, which means that Pes Axyx) presents, 
on +º, the same distribution as 2. 

Let us make q stationary in the Lagrangian way on 
tho set of the disc eta fields which satisfies equations 
(99) and (107) and the conditions 


Ult) = ul!) (115) 
Uintlo) = Usnlts) (116) 
Pull) = pnlti) (117) 
Palto) = pato) (118) 


We start by determining 5 T: 


dT=0)5 | vvpdA' = 025 IR Vibmp dA vi vim 
Es ae e A” 


— Ê, ; ) vob Ap dA, vim +52 ) v;ôv;p dA“ 
E A” E A 
(119) 


Considering (96), é T becomes 


OT = 2 U; ÔU, dA“ (| 20) 
ge 


Ncw as equation (49) is still valid, (120) can be trans- 
formed into 


OT = 5) v;ôv,;p dA* = 
e JAº 
=5L). v; Ou; p dA — EI. dus (v;p dAº) 
d 
= di 


v; du;p dA* — 
e 


2]. du;0;p dA* — :]. durão (po dA*) 


(121) 


As, by virtua of (90a), (90b) and (96), 


148 


:Q 
Db, Õ| o: dA“ =— 
a À Wy va A? dA“ Vôulim = = |) (122) 
we finally obtain 
Ez iz 
| ôTdt = [z] v;ÔU;p a8º| | XE du,0,p dA* 
ar e Y Fita 5 (123) 
or, by virtua of (115) and (116), 
ta ta 
| oT dt = -| Z| du; vp dA (124) 
Hi oe VA” 
As 
Pp. 
dy = ao it) (125) 
equation (124) vyields 


'2 iz RR 
[ ôTdt = -| z) pau) du, dA“ 
fi ty E el Ut 


a À (CrontPpolin + d;ponmi Prolint ja! OU di 


(|| 
where 
CrQLN = Z Probnti dA“ (127 
E A” 
d,poNML => ProW nba dA“ (128) 
e A” 


The expression of the internal energy variation can 
easily be determined if equations (97) and (98) are 
considered. We find 


=>) Wiwyôp dA'niy = 


e WA? 


= e | viyyp div ôdudA'r;y 
E A” 


== a PI GM = |. PY dA Ou; 


(129) 
or 
dE = e; ÔUIN (130) 
if we make 
ENS — > py dA“ (131) 
E A” 


Let us consider finally the terms of WV due to the 
external forces, 
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VK=—> | foxip dA + [ E) x,p dB (132) 
e “Ae B; 1P 
We assume /, and (p/p) to be of the form 
fi=Wnfin (133a) 
[r] |? (133b) 


D) = vs), 


within 4º and on Bm Bº 
Therefore, by virtue of (96), 


OV. = — | fiôdu,p dA + ba 0) du,p dB 
E A? e Bra B, P 
(134) 


so that, making 


A -5| usas — | [vim der] 
A” Beak, b! 


E 
we can write 


ôV. = —fin Ô4in (136) 


Consideration of (126), (130) and (136) vields 


E+ 
Op fm | | a C PON LÊ prin Eee d;poNM LÊ ro! NUM + 


=+- Ci +fa) du; di (137) 
The necessary and sufficiant condition for py beco- 
ming stationa:;y is thus 


pralCronttin + diponmP pQUintjn)— Gil -J.=0 
(138) 
The number of equations in (138) is equal to the 
number of unknowns V;n- 
The set of generalized equations (99), (107) and 
(138), together with the initial conditions (117) and 
(118) and the boundary conditions 


Pin = VAN) (139) 


are sufficient for the determination of the unknowns 
Vin Pim and Tim 

Unknowns u;y can be determined from (123) and 
tha initial conditions (115) and (116). 

We point out that the solution of the discrete 
equations (137), (99), (107) and (123), together with 
the initial and final conditions (115)—(118) and boun- 
dary conditions (138) make functional q stationary 
under the conditions (99), (107), (123), (138) and 
(115)— (118). 

It is important to observe that, as equations (99), 
(107) and (128) (or (28), (31) and (36a)) and their 
boundary and initial conditions form a self-contained 
system which does not involve the displacements, the 
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distribution of the velocities, density and pressure can 
ba dete:mined on A without the concept of displacement 
even being considered. 

Equations (129) are practically the same as the 
ones presented by Oden.” The continuity equation is 
different, however. On the other hand, Oden considers 
only one kind of nodes and not two kinds as we are 
obliged to do in the present paper. 


Order of approximation 


The elements of space À are now the inertial fields 
on 4. The norm of an arbitrary element xe X is given 
by 


(140) 


] 
|x|| = Jz (vv,;+p” +pPjdas | 


E A 


According to the minimum theorem, the subset 
Cc X is the set of all the continuous fields (i.e., 
tha set of all the fields which satisfy the continuity 
equation, the equation of state and the velocity boun- 
dary conditions) which satisfy the initial conditions 
(53-56). 

The elements of X' are the inertial fields which 
satisfy equations (90a), (91)— (93) within each subdo- 
main, and the elements of C' are the elements of X' 
which satisfy the generalized continuity equation (99), 
the generalized equation of state (107), the boundary 
condition (139) and the initial conditions (115)— (118). 
The B-projections of each element x of X in X' are the 
elements of X' which present the same nodal values as 
the displacements, velocities and density as x. 

The approximation theorem mentioned in section 
“Results of the theory of variational methods” states 
that a definition of distance can be introduced such 
that (1) holds, but the physical meaning of such dis- 
tance is not obvious in the general case. 

In what concerns the present case, a discussion 
of the physical meaning is intended to be made in the 
near future. Meanwhile, it does no harm to think of the 
distanco between two fields as the norm of the diffe- 
renco between them, i.e,, 


dix, x2) = x, —x1| (141) 
The variation of the functional y from an element 


xe X to another element within a very small distance 
from x is, considering (42) -—- (44), given by 


d,P = Ollô,ul)+ O(jó,v|) + O(d,p) 


d,P is thus of the order of ["*!if 

(i) nris the degrees of highest degree complete 
polynomial of arbitrary coefficients contained in the 
expressions for the velocities and densities within each 
subdomain Aº, 

(li) all the derivatives of the velocities and densi- 
ties associated with x. with order nl or less are 
bounded within each subdomain As. 

Equation (93) ensures that conditions (i) and (ii) 
are satisfied for n = 1. 


(142) 
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No reference is needed to the displacements in 
conditions (i) and (ii) because, as tha velocities remain 
the derivatives with respect to time of the displace- 
ments, these conditions are satisfied by the displace- 
ments whenever they ara satisfied by the velocities. 

Observing that  d(s,5,) is of higher order than 

v(ó,p). we can write (see Tu 
dº(s. 5) < |ô,pl+lô,9] (143) 
and state the distance between the exact and appro- 
ximate solutions is of the order of |!"* !/'2 (see section 
Viscous non-inertial flow') if: 

(a) condition (i) is satisfied, 

(b) condition (ii) is also satisfied by the exact 
and the approximate solutions (whatever the value of |), 

(tc) both the exact and approximate solutions are 
stable. 

Requirement (c) results from the fact that the 
variational principle is a minimum principla only in case 
of stability. 

We observe, however, that the exact and approxi 
mate solutions which are being considered refer to a 
problem which is not the usual one. 

The initial and final values of the displacements 
and density are indeed assumed prescribed (first pro- 
blem), while the values which are usually prescribed 
aro the initial values of the displacements, velocities 
and density (second problem). The first derivative of 
the density with respect to time is not mentioned be- 
cause its values is connected to tha values of the velo- 
cities and density through the continuity equation. 


In this second problem, condition (b) must hold 
both for the exact solution s and for the discrete solu- 
tion s, corresponding to the same initial and final 
values of the displacements and density as the exact 
solution, in order that it may be ensured that the. dis- 
tance between s and the approximate solution s, be of 
the order of 11" 1h= 

A proof of this statement can easily ba made by 
considering the Banach space of the generalized cons- 
traints &* and the linear function y'(x') which associates 
to each element of X' an element of 6.º As |s—s,|| 
is clearly of the order of/'"* ''2 and y' is assumed con- 
tinuous, it follows that [7 [B(s)]-v(si)le is also 
of the order of [Int 12 

On the other hand, B(s) and s' 
ned,, i.e., 


are isoconstrai- 


"LB(S)) = (8,5) (144) 


E 


Therefore, as operator y' is assumed bounded, 


7 (SD) — (solo: = OW" * 1/2) (145) 

Calling now "'p' ;s the subset of y with tha same 
range as y' but whose domain is the set Xp: Of the 
elements of X' which equilibrates the same system D' 
of external forces as s, , we can write, considering 


[| 
that s, belongs to X, and that y' is linear, 
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Iyp(sa— Sole = OU!) (146) 
Now, by virtua of tha uniqueness assumption, func- 
tion Yp' has an inverse and, by virtue of stability, such 


an inverse is continuous, so that (145) implies 
Isis! || = O(1M* 12) (147) 


and, therefore, 


|s=si;ll = 0(!!t 2) (148) 

Of course, it cannot be forgotien that the only 
error which is being discussed is the discretization error 
due to the use of the finite alement method for the 
integration of the continuous equations with resvect to 
the space co-ordinates, not with respect to time. 

The discreta equations are thus supposed to be 
exactly integrated with respect to time and no round- 
-Off error is considered. 


VISCOUS INERTIAL FOW 


We assume that a viscous inertial fluid is a com- 
plex body resulting from the association of two simpler 
bodies: a viscous non-inertial fluid and an inertial non- 
-viscous fluid. The magnitudes referring to the first 
will be denoted by index V. The magnitudes referring 
to the second will ba denoted by index /. The symbols 
with no index refer 1º the complex body, 

The kind of association can be defined by requi 
ring that the values of the stresses in the comolex body 
fjno index) result from adding the stresses of the sim- 
pler ones, and that the kinematical magnitudes (dis- 
placements, velocities, strains) and the density have 
the same values in the complex body and in the simpler 
ones. These kinematical magnitudes and the density 
are denoted therefore by no index. 

By virtua of the conventions just established equa- 
tions (8)—(11) may be rewritten in the form 


Cu Potal; =0 (8) 
o; = Quê, (10) 
p= —uê (11) 
and equations (28), (31) and (36a) in the form 
pij+p E pf, (28) 
RO 0 3 
= adivv = 31 
di f (51) 
p' = fp) (32) 
Adding equations (8') and (28'), we obtain 
Giji— (p' + pp), + pt! + [)- pó; = 0 (149) 


TECNICA 432 


Conside. ing that 


NE = | É 
O, é 0, (150) 
' | " 
Pp =p (151] 
f;+f; =, (152) 
(149) can be transformed into 
Gia Ps +ply— pb; =0 (155) 


Equations (153), (9'), (31), (10') and the equa- 
tion 


p=p)-xé (154) 
which results frem adding equations (11') and (32), 
ara tho continuous equations governing the viscous 
non-inertial flow. 

The 3 equations contained in (153), the 6 equa- 
tions contained in (9'), the 5 independent equations 
contained in (10'), equation (31') and equation (154) 
mako 16 equations. 

The unknowns a:e also 16: 3 velocities, 6 rates 
of strain é; and é, 5 independent stresses Gij. the pres- 
sure p and the density p. 

Eliminating the stresses and the rates of strain, 
wo obtain Navier-Stokes equation : 


l ' as 


Combining boundary conditions (14) and (15) 


with (37) and (38) ylelds 


õ 
Er —Én, = T, on B, (156) 
f [é 
p= on B, (157) 


Combination of the internal boundary conditions 
(16) and (17) with (39a) and (40) vields 


(0)'n—p'n; = (6;) n—p”n; (158) 
on T 


D; = UV; (159) 


Let us assume now tetrahedral subdomains of the 
same kind as the ones consideced in section 'Inertial 
ncn-viscous flow”. 

The discrete equations for the viscous non-inertial 
fluid area 


E .. | 
Mimnlim = Sn (23) 
Tho discrete equations for the inertial non-viscous 
flow are 
Post ind Di + dd. Tabu J=e! api = f) 
EroteronLtiN TE;ponmLProrint;m Til 
(138) 


* This means of course that the displacements, velocities and density take the same values in x 
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(99') 
(107') 


armpoPro + DimponProlin = O 
; | ro 
gimpoTim = hpolp) 


Cembining (23') and (138') yields 


E a % n al k | 
ProlCponLCin + diponmif'po! nt;m)— Cla = 4 


(160) 
where 
fu =Jutli= 2 | Wi fi dA + | Wsã, ar] 
| Lo Bea BH, 
(161) 


Equations (160), (99') and (107') form the system 
of discrete equations which governs the viscous inertial 
flew. lts solution provides of course an approximate 
solution for Navier-Stokes equation. 

Let us discuss now the order of approximation. 

A set of velocity, strain rate, density and pressure 
fields will bo called a viscous inertial field. 

The set of all possible viscous inertial fields forms 
a vector space which becomes a Banach space X if a 
norms is introduced. The norm of an element xe X 
is defined as 


v| = JzJ et ty? + pda] (162) 
E A” 


Let us consider now the Banach spaces X" and X' 
of all possible viscous and inertial fields and let xr 
and x denote the fields which correspond to x when 
tha decomposition of the comolex body into simpler 
ones is considered.* 

- We denote the norm of X”and X' by |x||, and 
x|ly. As 


|xly = Jz Eijês) an] (163) 
e A” 
|x|l = Jzf ento + ripar] (164) 
e “A” 
it is clear that 
x] < Ixiy+ xi; (165) 


We call the distance between two fields in X the 
norm of their difference, i.e., 


dlx, X2) = |x,— x5|| (166) 

The systems of external forces which can act on 
the body are defined by a body force density vector 
function f(x, ) and by a traction vector function m(X, [) 
on B,. 

We introduce the Banach space 4 whose elements 
are pairs of functions f(x. t) and f(x,t). The norm of 
an element De 4 may be defined as 


| x" and x. 
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pio= Ji[ xl] E fadas | K.xdB” | 
Pipor A Ho nB; 


(167) 


We also introduce the function 5(x) with domain 
X and range & which associates which each element 
of X an element of € such that equations (153) and 
(156) are saisfied. 

We assume the portion of the fluid which occupies 
the domain 4 decomposed into Eulerian finite elements 
and call X”., X! and X' the subspaces of X”, X! and 
X which correspond to such decomposition. 

Wa define the interpolation operators B” and 8! 
as the interpolation operators in sections “Viscous non- 
-inertial flow' and 'Inertial non-viscous flow. The inter: 
polation operator B is of course defined to be such that, 
if x” and x! represent the decomposition of x, B” (x) 
and 8! (x!) represent the decomposition of Blix). In 
other words, the B-projections in X' of the elements of 
X are the fields in X' which present the same nodal 
values of the displacements, velocities and density. 

We call: 

C” the set of fields in X 
crete boundary conditions (22), 

Cl" the set of fields in X which satisty (22) or 
(138') and also (99') and (107). 

C' the set of fields in X' which satisfy all these 
equations simultaneously, 

(4 the Banach space whose elemcnts are sets ha 


which satisfy the dis- 


functions Jilt) defined on the interval (to to); 
norm in £& being given by 
|D'lg: = JM futa. ) (165) 


de lx) the function with domain C' and range 4 
which associates to each element of C' the element! 
of “4 such that equation (160) is satisfied, 

ôL(x”) and delx! ) the corresponding functions 
for the viscous non-inertial and for the inertial non-vis- 
cous fluids, 

A”, A! and A the approximation operators corres- 
ponding to the two simple models and to the complex 
model which are being considered, which means that, 
if s” represents the solution of the viscous non-inertial 
continuous equations, A” (s”) represents the solution of 
tha corresponding discrete equations, and so on. 

As the type of association is such that the exter- 
nal forces acting on the complex body represent the 
sum of the forces acting on the simpler component 
bodies, we can write, if x” and x! represent the de- 
composition of x. 

rx) + dd!) = dA) (169) 

Assuming that the approximation conditions in 
sections “Viscous non-inertial flow” and 'Inertial non- 
-viscous flow' are satisfied by the expressions for the 
different fields within the subdomains 4º and by the 
pairs of solutions s” and A” (sv), s! and A! (s!], we 
can write 
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18 = AP (Sh = "OU E] (170) 
[88 — AS!) = OU 1) (171) 
and, Of course, 
CÁ PE RS id 172) 
= Bis, = O +=] (173) 


Combining (170) and (172), (171) and (173) vields 


(174) 
(175) 


1Br(S")— AM(s "ll = Oto 15) 
|B/(s!)— A!(s) ll = 04!" !"5) 


Now, as B“(s”) and B'(s!) represent the decom- 
position of B(s), we can write (see (170)) 


SL LB'(s")] + d4.[B'(s!)] = de[ Bis! 6) 
On the other hand, as 


| 


S"(s) + 05!) = (8) (177) 
there results, by virtue of equations (25), (135) and 
(161), 

SETAS") + LA! (SN] = de LA(S)] (178) 


Subtracting (178) from (176), and considering 
that functions op SM and dc are linear, we obtain 


SE. LB" (s")— A!(s")] + d-[B'(s”) — 


— A!(s9] = de[B(s)— A(s)) (179) 


But as such functions are also continuous, we can 
write, considering (174) and (175), 


6 [B"(s")— AM(S lg = O(”* 112) (180) 
WôLIB's!)— A!(s ll = O(!"* 12) (181) 

and then, considering (179), 
Ide [B(s)— Als): = OU" 15) (182) 


The uniqueness and stability assumptions enable 
us to conclude that 


| B(s)— A(s)|] = O(M* "* (183) 
Às, on the other hand, by virtue of (165), (172) 

and (173), 
lis— B(s)] = O(1!"* 15), (184) 

there results 
|s— A(s)] = O(1* 122) (185) 


o., considering (166) and the fact that Als) iss... 


ds. s4) = Ot t 2) (186) 
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The distance between the exact and approximate 
solution is thus of the order of |"! !)2 ijf the require- 
ments for convergence are met by each of the exact 
solutions s” and s!, which result from the decompo: 
sition of s, and by the discrete solutions A” (s”) and 
A! (s!). The latter, however, do not represent the decom- 
position of A(s), but the discrete solutions generated 
by the finita element method applied to the analysis 
of the viscous non-inertial and inertial non-viscous bo- 
dies acted by the same forces which act on the simple 
bodies within the frame of the continuous theory. As 
the distribuition of the external forces by the compo- 
nent bodies is altered by discretization, these forces are 
indeed different from the ones which act on the simple 
bodies within the frame of the discrete theory. 


CONCLUSIONS 


The application of a technique based on the decompo- 
sition of a complex body, whose mechanical behaviour 
is not governed by a minimum principle, into simpler 
ones for which extremum principles are known made 
it possible to apply the results of a previously develo- 
ped theory of variational methods to the discussion 
of the convergence of sequences of finite element 
approximate solutions of Navier-Stokes equation. 

The variational principles for the viscous non- 
-jnertial case are well-known. 

À variational principle for the inertial non-viscous 
flow, capabla of providing convenient discrete equations, 
is however one of the main contributions of this paper. 

The principle, which involves displacements as 
well as velocities, is Lagrangian, although the finite 
elements* are Eulerian. In other words, the decompo- 
sition of the domain into finite elements does not 
change with time, but the variation of the functional 
which the exact solution makes equal to zero is La 
grangian. 

A short discussion of the conditions under which 
such a variational principle is an extremum principle 
was included and a connexion was established between 
such conditions and the stability of the flow. 

Sufficient conditions were obtained for the dis- 
tance between the exact and approximate solutions 
being of the order of |" ' 1/2 | were | is the maximum 
diameter of the finite elements into which the fluid is 
assumed to be decomposed, 

Such conditions will be summarized now and 
practical convergencea tests will be searched, 

Assume that the expressions for the velocities, 
rates of strain and density within each subdomain à“ 
contain a complete polynomial of the nth degree. 

In order that the approximation error be of the 
order of [!"* !H2, mwo conditions must still be satisfied 
by the exact solution s and by the solutions of the 
viscous non-ine:tial and inertial non-viscous discrete 
equations cor-esponding to 


— the same forces which equilibrate the fields 
sY and s! resulting from the decomposition of s, 

— the same initial and final values of the displa- 
cements and density as s. 

These conditions are: 

(a) The derivatives of the velocities, rates of strain 
and density with order n + 1 or less are all bounded 
within each subdomain As. 

(b) The solution must be stable 

In order that convergence may be achieved, n must 
be ai least equal to zero. This means that the velocities 
rates of strain and density must be able to take arbi- 
trary constant values within each subdomain Aº, and 
tha: they must be bounded, as well as their first deri- 
vatives, within each subdomain, both in the exact solu- 
tion and in the discrete solutions mentioned above. 

Conce-ning the exact solution, physical intuition, 
a qualitative knowledge of the phenomena and general 
theo.ems on singularities can be used for establishing 
the boundedness of the different magnitudes. 

Concerning the discrete solutions, a supplementary 
sourco of unboundedness is connected with a possible 
heterogeneity of the types of nodes associated to some 
patterns of the mesh, as in the st.uctural case. 

This second source of unboundedness can be 
eliminated by resorting to patch tests similar to the 
ona which lrons* introduced in the structural case. 

A firs: test, devised to ensure the boundedness 
of the derivatives of the velocities end density, consists 
simply in checking for am arbitrary patch if the gene- 
ralized solution becomes exact when the boundary nodal 
values of the velocities and the initial nodal values of 
the velocities and density are the ones which corres- 
pond to an arbitrary permanent uniform distribution of 
the velocities and density and no external forces act 
on the body. 

That it is not necessary to make this test for each 
particular case is an immediate conclusion of the fact 
tha: equations (86) and (160) are indeed satisfied by 
uniform and permanent nodal values of the velocities 
and density and by fi = O. 

À second test, devised to ensure the boundedness 
of the derivatives of the rates of strain, is almost iden- 
tical to the one used in the structural case. It consists 
in checking for an arbitrary patch If the generalized 
solution becomes exact when the boundary nodal values 
of the velocities are the ones which correspond to an 
arbitrary linear distribution of the velocities, and the 
initial nodal values of the velocities and density also 
agree with such distribution and with a uniform ap. 

In the non-viscous case considered in section 
Inertial non- viscous flow” the boundedness of the rates 
of strain and their derivatives is not necessary for 
convergence. This second test is thus also not neces- 
sary. Attention must be called, however, to the facts 
tha: tho non-viscous flow is more easily unstable and 
that discontinuities of the velocities and density are 


* Only tetrahedral elements were considered with two kinds of nodes: vertices and mid-edge points. 
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possible without the rate of tha energy density beco- 
ming unbounded. 

It remains to be shown in a future paper how 
passing the patch tests is connected with the fact 
that certain magnitudes remain bounded when the di- 
mensions of the finite elements become indefinitely 
small. 
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Cálculo de estruturas reticuladas no espaço 


Analysis of space frames 


RESUMO 


Expõe-se sucintamente a aplicação do método de 
Cross na resolução de estruturas espaciais, constituídas 
por barras de momento de inércia constante ou variável, 
de eixo rectilineo, ortogonais entre si, ou não. 
Apresentam-se alguns exemplos de aplicação. 


Na nossa bibliografia técnica existem algumas pu- 
blicações sobre o tema que vamos abordar, designada- 
mente as indicadas nas referências [4], [5], [6]. [7] 
não havendo todavia nenhuma que conside:e o proble- 
ma com generalidade, 


A extensão do método de Cross, nos termos adian- 
to referidos, foi exposta por C, Catuneau, D, Anastesescu 
o |. Munteanu, na publicação referida em [1]. 


Como é sabido, no caso de estruturas a três di- 
mensões de barras perpendiculares entre si pode-se 
aplicar o método de Cross nas condições em que é 
utilizado cor-entemente para as estruturas planas, mas 
acrescentando a rigidez à torção. Neste caso, ao esta- 
belecer o equilibrio dos diversos nós, os momentos 
desiquilibrados são de flexão para algumas das barras 
a de torção para outras, havendo que considerar, para 
a determinação dos coeficientes de distribuição em 
cada nó, nas primeiras a rigidez de flexão e nas outras 
barras, a rigidez de torção, nas condições habituais. 
O mesmo não acontece com as estruturas de barras 
obliquas, e a dificuldade nestas resulta do facto da 
deformação de uma qualquer das barras que concorre 
num nó da estrutura, (por virtude da aplicação da um 
momento de flexão ou de torção no referido nó e que 
se distribua pelas barras que convergem nesse mesmo 
nó), ser simultâneamente de flexão e de torção. Não é 
possivel determinar, a priori, as componentes da rotação 
das barras que concorrem num nó, (condição necessá- 
ria para poder fazer a distribuição, pelas barras, dos 
momentos desiquilibrados), a não ser pelo estabeleci- 
mento da condição de compatibilidade inicial, que há 
que resolver à partida. 

Por esta dificuldade, mesmo a literatura estran- 
geira conhecida só tem, em regra, abordado, com o 
método de Cross, o caso das estruturas no espaço cons- 
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J. 5. BRAZÃO FARINHA 


SYNOPSIS 


li is succintly stated the application of the Cross 
method on the resolution of rigid space structures, cons- 
tituted by bars of constant or variable inertia, of recti- 
linear axle, orthogonal between them, or not. 

Somo application examples are presented. 


tituídas por barras perpendiculares entre si, onde tal 
dificuldade não existe. 


Para a resolução, no caso mais geral, de uma 
estiutura no espaço, há que adoptar dois sistemas de 
eixos ortogonais, um privativo de cada barra que con- 
corre no nó i conside;ado e o outro (constituído pelos 
eixos x, y, Z com origem no correspondente nó), de 
referência da estrutura. Do primeiro sistema de eixos 
dois deles (designados respectivamente por 1 e 2) coin- 
cidem com os eixos principais de inércia da secção e 
o outro (designado pelo algarismo 3) tem a direcção 
perpendicular ao plano da secção, conforme figura 1. 

É necessário orientar os eixos 1, 2, 3 das secções 
das barras que concorrem em cada nó da estrutura de 
modo que os respectivos sentidos coincidam com a 
convenção de sinais de Cross pa;ra os momentos das 
reacções exteriores (positivos quando dão rotação no 
sentido contrário aos dos ponteiros do relógio). O sen- 
tido a tomar para os momentos (vectores representados 
por setas rectilineas) será o da progressão da ponta 
do sacarrolhas quando se imprime a este uma rotação 
de sentido representado na figura por setas curvas. 
Porque esta representação é menos clara no desenho, 
que aparece em perspectiva, preferimos representar os 
momentos sob a forma vectorial. Uma convenção de 
sinais deste tipo, que não muda consoante se trata de 
extremidades esquerda ou direita das barras, tem a 
vantagem de definir claramente o sentido da actuação 
dos momentos instalados nas extremidades das barras, 
trate-se de momentos de flexão ou da momentos de 
torção. 

Pelo contrário, a orientação dos eixos x, y, z é indi- 
ferente, havendo apenas que garantir a coerência dos 
sinais dos cosenos e dos sinais das componentes dos 
momentos em relação a estes eixos. 
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Fig. 1 —- Sistemas de eixos de referencia e convenção de niveis 


positvos no nó | (Cross) 


Os ângulos que os eixos 1, 2, 3 relativos à barra 
i-j fazem com os eixos x, y, z são, respectivamente, 
designados por a, 5, y. São medidos entre sentidos 
positivos dos referidos eixos, e descritos do eixo x ou y 
ou z para o eixo 1 ou 2 ou 3. É a medida do ângulo 
assim definido, que permite determinar o valor do 
coseno do respectivo ângulo. 

Cada ângulo é referenciado pela indicada letra 
grega afectada de um índice constituído por duas letras 
separadas por um hifen definidoras da barra, seguidas 
de uma virgula de separação e de um número que desi- 
gna a qual dos eixos 1, 2 ou 3 o ângulo se refere, ver 
figura 1. 

Os momentos são definidos por vectores (os de 
flexão dirigidos segundo um ou outro dos eixos prin- 
cipais de inércia da secção, portanto 1 ou 2, em relação 
aos quais o momento de inércia | é respectivamente 
considerado; os momentos de torção são sempre diri- 
gidos segundo os eixos 3), e pelos sentidos da rotação 
no plano perpendicular ao vector. 
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| — Rigidez e coeficientes de transmissão 


Nas estruturas a três dimensões há a considerar 
as rigidezes junto das respectivas extremidades i, je 
segundo os eixos le 2 (de flexão) e 3 (de torção). 


f.1! — Vigas de momento de inércia constante 


No caso da barra de eixo rectilineo e secção 
constante perfeitamente encastrada à flexão e à torção, 
as rigidezes são iguais nas duas extremidades e os seus 
valores são dados pelas expressões 


em GUe E |, São os momentos de inércia da secção 
em relação aos eixos 1 e 2, respectivamente, e |. O 
momento de inércia à torção. 

Os coeficientes de transmissão tem os valores 
seguintes flexão n — 0,9 e de torção n, = -1. 

Os valores de |, |. e |, são indicados a segui:, para 
as secções mais correntes: 


Secção circular de raio r 


1º = 0,785 1º 


= 150 n=2,=2, 


Secção anelar de raios Fer TF; 


Secção rectangular de dimensões axb, colocada ao alto 
em relação ao eixo 1 (a <b, altura b) 


ab” ba! 


-— * ie 


' mi (| metes 
2 o q 


Secção rectangular de dimensões axb, colocada ao baixo 
em relação ao eixo 1 (a <| b, altu;a a) 


ba” ab” 
| —— dy E 
, 12 , 12 


| = RE 
| Ed | 1.630 Tr] + 0,052 [7] | | 
b b 
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- Go 


Rar 


Fig. 2 — Cálculo do valor de | 


numa secção metálica 


Jeso 
pa “Sa 


- «820,80, 


roj— 
A 


Fig. 3 — Barra de momento de inércia variável 


O quadro. 1 indica os valores de |, relativos à 
secção rectangular, para relações entre a altura e a base 
variaveis de 0,1 a 5. 

Para uma secção rectangular estreita, de altura b 


aºb 
muito maior do que a, |, = 


No caso de uma secção | ou T, tomar |, igual à 
soma dos valores que correspondem aos rectangulos 
em quo a secção se pode dividir. 

Assim, no caso da secção UNP14, fig. 2, para a 
alma da viga temos: 


0,73 x 14 
Dj E — 3 = 1,60 em! 
o para os banzos 
1X 
“4=2 fe = 3,53 


donde se tira: 
=", + 1" = 5,13 emt 


Nas secções de betão armado pode considerar-se 
apenas a área do betão, não entrando em consideração 
com a existência das armaduras de aço. 


1.2 — Vigas de momento de inércia variável 


No caso de secções de momento de inércia variável, 
vejamos, através de um exemplo, como se podem calcu- 
lar as rigidezes da viga representada na fig. 3, de 
largura constante a. 


Quadro 1 — Valores de 1,/1, 
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alt. da viga a | 
ae o Eq RA SS a 0,10 0,20 0,30 0,40 0,50 0,60 0,70 0,80 0,90 1,00 
larg. da viga b | 
1/1, 3,75 | 3,60 3,24 2,99 2,75 2,50 2,27 2,05 1,85 1,69 
It. da viga b 
Ato voa 2 | 420 | 140 | 1,60 | 180 | 200 | 220 | 240 | 260 | 280 | 3,00 
larg. da viga a 
FE 1,38 | 1,14 0,955 | 0,806 | 0,687 0,591 0,513 | 0,449 | 0,396 | 0,351 
alt. da viga b | 
si na iss 3,20 | 3,40 3,60 3,80 4,00 4,20 4,40 4,60 4,80 5,00 
larg. da viga a 
1/1, 0,314 | 0,282 | 0,255 | 0,231 0,211 0,193 0,177 | 0,163 | 0,151 | 0,140 
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A aplicação de um momento flector K, iq na 
extremidade h de modo que faça rodar o apoio h de 
rotação 6, = 1 conduz à instalação do reacções de 
apoio R e momentos M ao longo da viga de tal modo 
que se verifica o sistema de equações (*) 


Í ri dx al x 
a * 


| K pill | 


| O o 


sistema este que determina os dois valores 


Fa xe 
Í E | | dx 
Eis * 
e o ÃO, al es 
Codx [ x x E 
| | ox | -— F dx 
“o É Ch do a O ia 
al 
x 
| dx 
É la 
= Ee 
los 
| dx 


Estas expressões são válidas para a determinação 
de Kpjo O UNR trocando a altura da viga com 
a sua largura, Para o caso anterior a altura da viga é 
variável tb, ) ea largura (a) é constante, agora será 
a altura constante a e a la:gura variável (b.) 

Tomando para referência a altura b na zona de 
secção constante, 


calculando-se através desta expressão os valores que 
constam do quadro 2. Tomando depois, para referência 


a largura “a” na zona de secção constante, obtem-se 


o momento da inércia em relação ao eixo 1, na zona 
em que a viga tem secção constante. 

Péla regra de Simpson teremos, em relação aos 
valorea do quadro 2, 


>< [ 0,20354+ 1,00000+2(0,26340+0,34924+0,47683+ 
+0,67497 + 5x1,00000)+4(0,23090 + 0,30230 
+0,40644+0,564447+0,81629x5xX1,00000 ] — 
0,73357 
| 


to 


e, analogamente, 


x |- 
| dx = AR a é E 


Ú Ú 


x? E 
= dx = 0,31923 — 


E 


0,05 | 


Fa 


3 o 


E 
O 


> | 1,63399+42,77778+2 (1,78063+ 1,95618 2,17014 |- 
-H 2,43665+5xX2,77778)+4(1,70416+1,86428 + 
+2,05761+2,29568+2,59605 , 5X2,77778 |= 


I 
= 244172 —— 


x |z xe had 
7 — 1.32805 —: | = x = 0,91825 


lo * “em o 


Teremos portanto, da acordo com os valores do 
quadro 2, 


0,31923 |º EL, 


K41 2 >>> >>> 
h-i, 1 0,73357 1X0,31923 |'—(0,44853 |): 
El, 
=2067 —— 
I 
0,44853 
a=(->—— -— 1)= 0,404 
"hi, fee ) 


o dos valores do quadro 3, 


0,91825 | EI, 
Ko 2——— ————— 
N-i,2 244172 1x0,91825-(1,32805 |2): 


El, 


= 2,05 


jo= | — 1)= 0,446 
Mri,2 0,91825 


(*) — J. S. Brazão Farinha “Cálculo de Estruturas Hiperestáticas” Técnica n.º 23, Outubro 1953, pág. 15 


158 


TÉCNICA 432 


Quadro 2 — Cálculo da rigidez de flexão e do coeficiente de transmissão K, .,; wy.;j da viga de momento de 


inércia variável 


| Alturas Mom. de Relações 
E | Anadasás da viga inércia E 
Secções | 7 b. | , Ber 1 x xº 
b el b ) Es o | 7 
| E 

h 0,00 1,70 4,91300 | 0,20354 —— | 0,00000 —— | 0,00000 —— 
1 c,05 1,63 4,33075 0,23090 0,01155 | 0,00058 
2 0,10 1,56 3,79642 0,26340 0,02634 0,00263 
3 0,15 1,49 3,30795 0,30230 0,04535 0,00680 
4 0,20 1,42 2,86329 0,34924 0,06985 * 0,01397 
5 0,25 1,35 2,46038 0,40644 0,10161 0,02540 
6 0,30 1,28 2,09715 0,47683 0,1432305 0,04292 
7 0,35 1,21 1,77156 0,56447 0,19756 | 0,06915 
8 | 0,40 1,14 1,48154 0,67497 0,26999 0,10800 
9 0,45 1,07 1,22504 0,81629 0,36733 0,16530 
10 0,50 1,00 1,00000 1,00000 0,50000 0,25000 
11 0,55 1,00 1,00000 1,00000 0,55000 0,30250 
12 0,60 1,00 1,00000 1,00000 0,60000 0,36000 
13 0,65 1,00 1,00000 1,00000 0,65000 0,42250 
14 0,70 1,00 1,00000 1,00000 0,70000 7,49000 
15 0,75 1,00 1,00000 1.00000 0,75000 0,56250 
18 0,80 1,00 1,00000 1,00000 0,80000 0,64000 
17 0,85 1,00 1,00000 | 1,00000 0,85000 0,72250 
18 | 0,90 1,00 1,00000 | 1,00000 0,90000 0,81000 
19 | 0,95 1,00 1,00000 1,00000 0,95000 0,90250 
, | 1,00 1,00 1,00000 1,00000 1,00000 


1,00C00 


. É Ei nc 


Quadro 3 — Cálculo da rigedez de flexão e do coeficiente de transmissão K, 5: Myjo: 


da viga de momento 


do inéccia variável, para a/b = 0,60 


Altura Momentos ce Quocientes 
da viga iné.cia Ep 
Secções Abcissas D. | E ) D. 1 x x 
b “".b b ho | li; mA 
1 | à 
h | 0,00 1,70 0,6120 1, 1,63399 0,00000 —— | 0,00000 — 
| o o E 
| | 0,05 1.63 0,5868 1,70416 0,08521 0,00426 
2 | 0,10 1,56 0,5616 1,78063 0,17805 0,01781 
3 | 0,15 1,49 0,5364 1.86428 0,27964 0,04195 
e | 02 1,42 0,5112 1.95618 0,39124 0,07825 
E | 0,25 1,35 0,4860 2,05761 0,51440 0,12860 
6 | 0,30 1,28 0,4608 2,17014 0,65104 0,19531 
7 0,35 1,21 0,4356 2,29568 0,80349 0,28122 
8 | 0,40 1,14 0,4104 2,43665 0,97466 0,38986 
9 | 0,45 1,07 0,3852 2,59605 1,16822 0,52570 
10 | 0,50 1,00 0,3600 2,77778 1,38889 0,69444 
11 0,55 1,00 0,3600 2,77778 1,52778 0,84028 
12 0,60 1,00 0,3600 2,77778 1,66667 1,00000 
13 0,65 1,00 0,3600 2,77778 1,80556 1,17361 
14 0,70 1,00 0,3600 2,77778 1,94444 1,36111 
15 0,75 1,00 0,3600 2,77778 2,08333 1,56250 
16 0,80 1,00 0,3600 2,77778 2,22222 1,77778 
17 0,85 1,00 0,3600 2,77778 2,36111 2,00694 
18 0,90 1,00 0,3600 2,77778 2,50000 2,25000 
19 0,95 | 1,00 0,3600 2,77778 2,63889 2,50694 
; 1,00 | 1,00 0,3600 2,77778 2,717778 2,77178 


Notemos que estes valo;es não dependem da re- 
El 


lação a por ser uma constante que não afecta o valor 


das fracções que dão K e q. 

Da modo análogo ao anteriormente indicado, a 
aplicação do momento flector K-n na extremidade 
|, de modo que faça rodar o apoio i da rotação 6, = 1 
conduz ao sistema de equações 


dx |—x 
K, h,l | + R, hi | dx = E 
A ç ”, I, 
ij t 
É pg “=: 
Kona | dx + Mind | tr 8 
bi de I, 
O [a] 
donde se deduz 
| (1=)? 
| e dx 
“q E 
Ko ET De e 
= al dx ea)a DA RR ? 
— —d | | dx) 
) P . , E , 
t 0] ts 
1| 
[=x 
| | | a 
, 
] epa seir À 
“ih, q 41 
|I—x)= 
| ( ) dá 
“ , 
LR] 
Teremos pois, de acordo com o quadro 2 
1 = 
(I—-x)> 1 x: 
| Ermo dx = | —. UR O dx-2 1X 
a = no Ri = , 
MR] it LR] 
a) 
X =” [+ 
X | —— x — 0,15574 — 
ni X Ho 
it 
o! = a] « : 
| a dx = | e | dx = 
a , is ” pd d, 
0 LJ 13 
= 0,28504 — 
ha 
e portanto 
0,15574 El, 
K. | >———— = 
ih 0,73357x0,15574-0,28504º | 
El 
=471— 
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0,28504 


Ni-ni= 015574 —- | = 0,830 


Da modo análogo, será, na conformidade com o 
quadro 3, 


ml 
|—x)* 
| ra dx = (2,44172+0,91825-2X 1,32805) Y 
, |, 
[8] 
[3 E 
j-—— = 0,70987 —— 
ti IR) 
a Ie e 
| É dx=(2,44172-1,32805) = = 1,11367 —— 
. h, “0 ” 
[6] 
0,707387 EI El 
Kg DL TE É = 1a — 
=? 244172X0,70387-1,113672 | | 
1,11367 
—4 =0,582 


“i-h2  0,70387 


A aplicação do momento de to:ção K, E à extre- 
midado h da ba:ra h—i, instala um momento constante 
ao longo do vão, (donde se conclue que o coeficiente 
de t;ansmissão continua a ser—1 nas vigas de momento 
de inércia variável), e significa também que a rigidez 
do to ção da viga é dada por 


=+k| 


Ls 


a q? a 
|, sendo k = 4 | =) 1-0,630[-1-) + 
UR: b 


k 


sons) 


r | 


| 


portanto dependente das dimensões a, b da viga. Na 
hipotese de ser a = 30 cm, b= 50 cm (quadro 4), 
conforme fig. 3, será: 


0,05 | 7 
—— 0,52510 + 1,10925 + 2 (0,58717 + 


3 


| 0,66581 + 0,76861 + 0,90849 + 5 X 1,10925) + 
+ 4 (0,55441 + C,62404 + 0,71355 + 0,83278 + 


1º 0,051 
+ 0,99907 + 5 X 1,10925) | — —. Y 
| 3 
| | 
X 55,66741 — = 0,92779 — 
o Ho 
1,078 GI, 


“a=15> Ki-ns= 
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Quadro 4 — Cálculo da rigidez de torção K, . I, 


; da viga de momento de inércia variável, para 


a= 30 cm; b= 50 cm 


Altura 
Secções id da viga | 
b (em) 

h 0,00 | 85,0 
1 0,05 81,5 
2 0,10 78,0 
3 0,15 74,5 
4 0,20 71,0 
5 0,25 67,5 
6 0,30 64,0 
7 0,35 60,5 
8 0,40 57,0 
9 0,45 53,5 
10 0,50 50,0 
11 0,55 50,0 
12 0,60 50,0 
13 0,65 50,0 
14 0,70 50,0 
15 0,75 50,0 
16 0,80 50,0 
17 0,85 50,0 
18 0,90 50,0 
19 0,95 50,0 
1,00 | 50,0 


No caso da viga da altura constante 50 cm em 
todo o seu comprimento, teriamos um valo; pouco di- 
ferente do anteiior 


0,902 Gi, 


= &-n5"— 


Kn = 13 


1.3 — Módulo de elasticidade G 


O módulo de elasticidade transversal G pode er- 
primir-so em função do módulo de elasticidade longitu- 
dinal E pelas relações 
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| 

F EM 

E 
1,90437 1, 0,52510 7 
1,80370 0,55441 
1,70306 0,58717 
1,60246 0,62404 
1,50191 0.66581 
1,40143 0,71355 
1,30104 0,76861 
1,20079 | 0,83278 
1,10072 | 0,90849 
1,00092 0,99907 
0,90150 1, 1.10925— 

o 
0,90150 1,10925 
0,90150 1,10925 
0,90150 1,10925 
0,90150 1,10925 
0,90150 1,1C925 
0,90150 1,10925 
0,90150 1,10925 
0,90150 1,10925 
0.90150 1,10925 
0,90150 1, 110925 — 
o 
nÊ E 
> 2(n+1) aqi 


1 
em que n é o módulo de Poisson e r=——— o coefi- 
n 


ciente de Poisson. 

Da acordo com o anexo | do REBA, para o betão 
armado » = 0,20 donds E = 2,40 G. Para a construção 
metálica toma-se, em regra, »r = 0,25, donde E=2,50 6. 


1.4 — Barras articuladas num dos apoios 


Indicamos atrás os valores das rigidezes das barras 
de eixo rectilineo e momento de inércia constante en 
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castrada nos apoios. Há todavia a referir que, no caso 
das bar.as encast:adas na extremidade i e articuladas à 
to:ção e à flexão, respectivamente, no outio extremo 
|. as regidezes são exp-essas por 


3 EL, 3 EI, 


Nesto caso os coeficientes de transmissão da ex- 
tremidado i da barra paia a extremidade j são nulos, 
E dp — O. 

Notemos que uma barra pode, independentemente, 
ser articulada numa das extremidades à flexão (ou to:- 
ção) e ser nela encastrada à torção (ou flexão), haven- 
do que proceder em conformidade. 

Com efeito, uma viga rígida pode ser encastrada 
à flexão (mas não à to:ção) por uma parede esbelta 
disposta no seu plano, conforme fig. 4a). Analogamen- 
to em relação à fig. 4 b) quanto à possibilidade de en- 
castramento à torção, e não à flexão. 


| aii 


Fig. 4 — Encastramento em i, à flexão e à torção 
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1.5 —. Casos de simetria 


Há ainda que considera; os casos de possivel 
simplificação dos cálculos, por razões de simetria, uma 
ver que então se pode considerar uma parte apenas 
(metado ou um quarto da estrutura) como se o restante 
não existisse. Para as barras cortadas tomam-se então 
iigideze: virtuais K 

No caso de uma estrutura simétrica cujo plano 
da simetria corta pelo meio uma ou mais barras, de 
tal modo que a estrutura se comporta como perfeita- 
mente encastrada (ou como articulada) na secção média, 
tomam-se nas barras cortadas rigidezes virtuais à flexão 
K, =0,5 Kf e K, = 1,5 Kf respectivamente no 
caso de actuação de solicitações simétricas ou anti- 
-simét.icas, ver fig. 5. Para as rigidezes de torção tere- 
mos, também respectivamente, K, = 0;K, =2kKt 

No caso do plano de simetria da estrutura passar 
po: um ou mais nós, de tal modo que as bar.as que 
coincidem com o referido plano não vodem te: rotação, 
consideram-se perfeitamente encastradas as ba:ras cujas 
ext emidades se inserem no plano de simetria, 

Pelo contrário, nas barras que se situam no plano 
de simetria e que podem rodar, e são divididas em 
duas metades longitudinais, toma-se em cada metade 
das barras, metade das ca.gas actuantes, e uma rigidez 
virtual de cada uma das partes igual a metade da real. 
Isto porque a cada uma das metades das bar;as corta- 
das correspondem metade dso esforços totais nela insta- 
lados. No final, portanto, terá que se multiplicar por 2 
todos os valores encastradas para as meias barras. 


2 — Fundamentos do método 


A aplicação do método de Cross a uma estrutura 
no espaço parte da situação de encast:amento perfeito, 
isto é, da aplicação, em todos os nós, dos momentos 
do encastramento perfeito M que respeitam as cargas 
aplicadas nas diversas bar;as, 

A libertação de um qualquer nó “ij”, corresponde 
à aplicação de um momento de sentido contrário ao 
aplicado nesse nó, isto é, - Mº , uma vez que em 
cada nó equilibrado, EM = &L. O momento - = Mº 
distribui-so pelas diversas barras que concorrem no nó 
|, de uma maneira que depende do valor das compo- 
nentes 9 ,,: 8,7: 9,; da rotação do referido nó. Os 
valo:es destas são dados pelo “sistema característico” 
de cada nó da estrutura, a seguir indicado para o nó 
|, é no qual os somatórios se estendem a todas as 
diversas barras que concorrem no referido nó, 


e 


O 4 *kK cos" a+g vi = K cos acos dj + 


- 9,4; * K cos a cos y + Mep = 0 


É x, 


UR K Z cos 8 cos y + M,i =0 


| EK cos a cos 4 + 0447 K cos? 5 + 
1) 


9,42 K cos a cos y8,; + cos 8 cos y 


+ 0,,2K cos" y+M,, =0 
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a) O plano de simetria corta as barras transver- 
salmente pelo meio 


—- 


SOLICITAÇÕES 


ANTI-SIMETEICA SIMETRICA 


Kyç=05 kF 
kyt = O 


kvg =1,5 kg 


ACTUAÇÃO DE 


kyvt= E ky 


b) O plano de simetria corta as barras longitudi- 
nalmente pelo meio 


ARO 

10 « 
ue a» ( 

< S(kyg 05 kg 
0 

E = kyt= 05k4 
u Ú 

D ulkyp=askg 
q cu 

35 

e a kyt=05 ks 
- E á 

U = 

«< 


Fig. B — Simplificações por vritude da simetria da estrutu- 
ras (barras de secção constante) 


Nestes sistemas de equações, K, nrKi-nz Ka is 
representam as rigidezes, referidas aos eixos 1 2,3 
da barra i-h, uma das barras que concorre em i. Além 
desta barra, há a considerar todas as outras i-j; i-k etc. 
a que pretence o nó i. 


Analogamente, Cp Cn? Tn CIC. represen- 
tam os ângulos que os três eixos da secção da barra 
i-h fazem com o eixo dos x e de modo análogo Bina 
Bi no Bio h,3 (com o eixo dos j) e 7. np: Vinha, 
Y1-n3' (com o eixo dos z). Além desta barra há que 
considerar todas as outras barras a que pertence o nó i. 


Finalmente, Mir Ma M, são as compo- 

nentes da resultante dos momentos aplicados ao nó i, 
segundo os eixos x, y, z, entrando com todos os mo- 
mentos aplicados a cada uma das barras que conco- 


rem nesse nó. 


Conhecidos os valores de 4! para o nó i) sistema 


de equações 1) as componentes M' dos momentos 
desiquilibrados (numa 1.º distribuição) e que se aplicam 
na extremidade i da barra ih, serão obtidos através 
das expressões 


Mini Kona | cos aj nata À 


x C0s Bj na 9!'zi COS fina | 


M'i-n2 = K- na) é xi 008 q n2 +60,4X 
2) 
x vos Bj n2 7 8º,4005 tina | 


M'i-h3 — S-na) 014; C08 4, n3 TO, À 


XcosB;. na T 8,00 Yi. a | 


Destas componentes, transmitem-se para a outra 
extremidade h da barra i-h os valores 


Ti— ha Mn * TMi- ha M,. h,2 M', -n3 


Do mesmo modo, para as restantes barras que 
concorrem no nó i. O indicado repete-se com valores 
92, 6% etc. até que as correcções introduzidas sejam 
consideradas desprezaveis, como é correntemente prati- 
cado nas estruturas planas. 


Em cada libertação de um nó |, há, pois, que 
resolver um sistema de equação (1) e depois temos de 
aplicar as expressões (2) para cada uma das barras 
que conco:rem no referido nó i, afim de determinar as 
componentes dos momentos desiquilibrados (segundo 
os eixos 1, 2, 3, de cada barra). Depois, em todas as 
barras, uma transmissão de momentos para a extremi- 
dade oposta a |. 

No caso de uma estrutura de barras ortogonais 
entra si, escolhendo um sistema x, y, Zz de referência 
do nó em coincidência com os eixos das duas ou três 
barras que nele concorrem, é sempre ZK cos a cos 8 = 
= EK cosa cos y= EK cos 8 cos y = O donde resul- 
ta que o sistema (1) se transforma no sistema (3) com 
a solução (3') imediata. 


0,17 K cos? a + M,|= O 
3) 9,/ZK cos? B+M,= 0 


A T cos? y + M; | = (6) 


(+) Se, neste caso, não actuarem forças nas barras cortadas longitudinalmente, consideram-se as barras perpendiculares a 


estas encastradas nas extremidades, e como se aquelas não existissem. Ver exemplo 4.º para esclarecimento. 
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[E] RR M e 
x, | = K cos" « 
M.. 
Vo 
3) oe ———— En 
Vl = K cos” B 
em M,, 
o” ES = a 
| ' z K cos* y 


Notando que para o caso de barras perpendiculares 
o valor de cos? é sempre igual a 1 (ou a O) deduz-se 
imediatamente, por substituição destes valores do siste- 
ma 2), que para este tipo de estruturas é generalizável o 
método de Cross nos termos em que é correntemente 
aplicado nos sistemas planos. Só temos, como já atraz 
dissemos, que considerar para cada barra a rigidez de 
flexão ou de torção, consoante o momento desiquilibra- 
do é de flexão (e portanto de flexão ou torção para a 
barra perpendicular, conforme se trata de barras compla- 
nas ou existentes em planos perpendiculares entre si). 
Ou então o momento desiquilibrado é de torção para 
uma dada barra, e de flexão para as que lhe são per- 
pendiculares. 


3 — Exemplos de aplicação no caso de estruturas 
de barras perpendiculares entre si 


Exemplo 1.º — Determinar os esforços actuantes 
na estrutura constituída por duas vigas simetricamente 
cruzadas, de betão armado e secção constante, e que 
apoiam uma laje na qual actua a carga uniformemente 
distribuida g (kg/mº), conforme fig. 6. 


A resolução do problema é considerada em duas 
etapas, na 1.º flexão nas barras h-i; i-h'; j-j; i-j/ com 
o nó central | fixo, quanto a deslocamentos na vertical. 
Do acordo cm a fig. 7 o nó i está em equilibrio, (mo- 
mentos M, iguais nos dois sentidos) não actuando 
portanto momentos de torção, dada a simetria da estru- 
tura. 

Na 2.º etapa considera-se o deslocamento A, na 
vertical, do nó i, que dá lugar a momentos M", também 
equilibrados. 


Não há portanto que fazer qualquer distribuição 
de momentos. 

O valor de A (flecha elastica em i) determina-se 
pela condição =V = O de equilibrio do nó i. 


1 | Bos 33 
1 A | PN em 
à Pa Ta? 184 P 
24E| 24EI 933 - 
Te = + ca E à = >= q 
E E 184 
11 pl! 
— 4472EI 


Os momentos flectores actuantes em i serão 
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33 25 
o 368 92 


| pl: gls 
fra ia 
] 16º 32 


nos dois sentidos, isto é, nas duas vigas que em i 
se cruzam. 


Fig. 6 — Estrutura simétrica, constituída por duas vigas cruzadas 


de igual vão, que sustentam uma laje de betão armado 


Exemplo 2.º — Determinar os momentos no pórtico 
indicado na fig. 8, no caso de actuar, normalmente ao 
plano da estrutura, uma carga uniformemente distribuida 
(vento) de 200 kg/m2, Considerar E= 2.00 G e as 
secções transversais das barras, rectangulares de lar: 
gura constante 0,30 m. 

Teremos nos montantes g = 200 X 0,5 = 100 
kg/m e na travessa g' = 200 X 0,72 = 144 kg/m. 

O problema é resolvido em duas etapas, conside- 
rando na primeira a rotação dos nó, sem deslocamentos 
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Fig. 7 —- Representação dos esforços que actuam no nó | da 
estrutura do 1.º exemplo 


no sentido perpendicular ao plano da estrutura, s consi- 
derando na segunda, etapa o deslocamento A dos nós 
i e j, no sentido perpendicular ao plano do pórtico. 
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Os momentos de encastramento perfeito corres- 


pondentes serão 


nos montantes 


100 x 5,02 

Mº = RR = 208,3 kgm para a carga 
GEI 

Mº = Eq a = 0,24 El A para o deslocamento 


na travessa 


144 XxX 10,0º 


| Mº = —————— = 1.200 kgm para a carga 


12 


| Mº= o para o deslocamento A, por razão de simetria 


Po: esta mesma simetria tomamos metade da estru- 


tura o adoptamos para rigidezes da barra i-ij dl, 


T2 XxX 30º 30 XxX 723 
=> = 182.000 em3; |, =———=" = 
12 12 

= 933.120 cmg), 
rigidez de flexão 
2 E 
| «=X 162.000 = 324 E 
Si-ia 1.000 
2 E 


K. —: 
i-12 1.000 


rigidez do torção K| ,=0 


30 X 50º 
| th E = = 
Para a barra i-h (1, 12 

50 X 30º 
= 112.500 cmd; |, = Sn 


= 312.500 cmd; |, = 2,50 |, = 


= 281.250 cm4), será 


Fig. 8 — Estrutura do 2.º exemple 
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yi 
Õ = — =— FE e —— 
y w K cos? A 900 
M. 
Bam — — Es. = 0 


W K cos* Y 


Os momentos da distribuição a aplicar na barra 
h—i são, no vértice | (quadro 5) 


Quadro 5 — Valores característicos do nó i (2.º exem- 
plo). Eixos de referencia da fig. 9. 


Fg. 9 — Sistemas de eixos de referencia e distribuição dos 
momentos correspondente ao 2.º exemplo 


rigidez de flexão 


4 E 
Kn1= 600 * 112.500 = 900 E 
4 E 
K-na= so X 312.500 =2.500 E 
Gi, 
igidez de torção K, O —— 2 
rigidez Ç i-h3 = 500 
2,5 X 112.500 E 
=———— = 225 E 
2,5 x 500 
Escolhemos os eixos de referência x, y em coin- ML, 900 x ——— = 208,3 kgm 
cidência com os eixos das barras da estrutura, e o eixo 900 
erpendicular ao plano desta, fig. 9. 
paia E E ; 1.200 
M. .=-— 225 X = — 490 kom 
Para a 1.º etapa teremos no nó i aplicados os i=h,5 549 ” 


momentos hM, = — 1200 (*), M, = — 208,3: M, —o. 


o no vértice h os momentos transmitidos, ou se) s 
De acordo com 3' será t ou seja este 


valores multiplicados pelos respectivos coeficientes de 
transmissão, 0,5 no 1.º caso por se tratar de flexão 


Mi 1.200 
E eMi=— Sa O 4. leixo 1) e — 4,0 no 2º caso por se tratar de torção 
ad a = K cos: EF 549 (eixo 3). 
(Continua) 
(*) Notemos que dada a orientação dos eixos da referencia, este momento pode exprimir-se por duas maneiras, M, = —1200 


ou Mi; |= + 1200 hgm 
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Limite superior do número de Froude em canais 


prismáticos e escoamento uniforme 


SUMÁRIO 


O auto: estuda a variação do número de Froude 
em canais prismáticos e escoamento uniforme e a possi- 
bilidade de o mesmo admitir um limite superior. Apre- 
sentam-se, igualmente, algumas aplicações numéricas. 


1 — INTRODUÇÃO 


O número de Froude é um parâmetro importante na 
caracterização dos escoamentos em superfície livre de 
fluidos incompressiveis. A variação do referido parâme- 
tro está dependente do andamento das curvas ds declive 
crítico em função da altura crítica. Em determinadas 
circunstâncias, o número de Froude apresenta um limite 
superior cuja expressão se deduz no presente estudo O 
qual se refere a canais prismáticos de leito fixo e a 
fluidos incompressíveis monofásicos, em escoamento 
uniforme. 


2 — NUMERO DE FROUDE EM ESCOAMENTO UNI- 
FORME 


2.1 — Secção rectangular 


No caso de escoamento em canais com secção 
rectangular, o número de Froude é definido [1] pela 
seguinto expressão 


U 

F = e 
/ A 

V mi, (1) 


a qual, se transforma, no caso muito corrente de se 


pode: considerarA = a = 1, na conhecida expressão 
U 
F = ES 
Voe h, (2) 


No escoamento uniforme, a altura poderá ser igual, 
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ANTÔNIO BETAMIO DE ALMEIDA 


SUMMARY 


The author studies Froude number range in pris- 
matic channels and uniform flow and the possibility of 
an upper bound for it. Numerical aplications are also 
shown. 


menor ou maior do que a altura crítica, consoante o de- 
elivo do canal for, respectivamente, igual, superior ou 
inferior ao declive crítico, Num canal com determinado 
declive, largura e rugosidade o regime de escoamento 
uniforme poderá depender do caudal escoado e da fór- 
mula de resistência utilizada. No caso do escoamento 
ser puramente turbulento, as leis de resistência mais 
frequentemente utilizadas [2] [3] são as que se indi- 
cam no Quadro 1. 


QUADRO 1 

LEI DESIGNAÇÃO EQUAÇÃO 

1 Chézy V=C, VRI 

2 Chézy U=C VRi 

3 | Gauckler - Manning U=KR25 il? 

| | 14.8R 
4 Karman - Prandtl |U = 5.6 log 4 V gRi 
e 


“À aplicação das leis apresentadas conduz às expres- 
sões para o declive crítico que se mostram no Quadro 
2. O declive crítico é função do caudal por intermédio 
da altura crítica, através da seguinte expressão que 
minimiza a correspondente à da energia específica : 


BP g (3) 
Sempre que, para um determinado caudal, o valor 


Sa. inferior a |, o escoamento uniforme é rápido e 


o correspondente número de Froude superior à unidade. 
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Na figura 1 apresentam-se as configurações da função 
il. NOS casos 1, 2 e 3. Nas figuras 2 e 3, relativas aos 
casos 4, 5 e 6, os valores de e são apresentados em 
função de ho e do parâmetro B, para um canal com 
secção rectangular para o qual se utilizou o valor 
K = 70 a que corresponde o valor de e = 2.6 mm for- 
necido pela expressão 


26 16 
e=( ç] (4) 


citada em [2]. À lei 2 comporta-se dum modo análogo 
ao da lei 3 sendo, contudo, a equivalência entre os 
parâmetros da rugosidade dependentes da altura h,.- 


O traçado das curvas foi efectuado automatica- 
mente mediante um programa escrito em linguagem BA- 
SIC para um computador H. P. 9830 - A, munido de 
«plotter». 


Das figuras podem-se extrair algumas considera- 
ções importantes, a saber: 


QUADRO 2 
Caso Lei Secção rectangular Declive crítico Figura 
1 1 Largura infinita 9 4 
(caso de Bresse) Cc, 
ho 
E 1+2— 1 
2 1 Largura finita B Cc; | B ) 
lá 
ue g 1 
3 3 Largura infinita ge [— 1 
Kº 4h. 
1 [ê 
h 415 1 15 
4 3 Largura finita B —— 1+2— 2 
| E B h 
| L 
1 á h 
5 4 Largura: infinita 32 log | 14.8 E 3 
| 9/ 14.8 ho 
ini | en creio | 
a 4 Largura finita B 32 g & 452 ck ) 3 
B 


NOTA — O valor 14.8, nos casos 5 e 6, foi obtida sem terem conta a correcção proposta por De Marchi [2]. 


le 
Caso Ê 
4 
A > 
ã; | = Coso! 
e | E 
Ecso3 
he 
Figura 1 — Andamento das curvas | =i ” (hc) nos casos 1, 


2 aegs. 


— No caso 1, o tipo do regime uniforme fica defi- 
nido independentemente do caudal escoado. 


— No caso 2, o tipo de regime varia consoante o 
caudal. Verifica-se a existência de um declive minimo 
(g/c;) abaixo do qual o regime é sempre lento. Para 
declives superiores o regime começa por ser rápido, 
tornando-se lento em virtude de aumento do caudal. 


— Nos casos 4 e 6 verifica-se a existência de um 

valor mínimo em do declive crítico para um determina- 
do intervalo de B. Num dado canal (secção, rugosidade 
a declive constante), o regime do escoamento uniforme 
vai variando consoante o caudal (figura 4). 
Para O < Qio regime é lento (F < 1), para OQ entre 
QieQso regime é rápido (F > 1), para O > Os o 
regime é lento (F< 1) epara O = Qio 0=05s0 
regime é crítico (F = 1). 


— Quando a largura tende para infinito, tende-se 
para os casos 3 e 5, verificando-se que O s tende, igual- 
mente, para infinito (figura 4), 


— Existem diferenças entre os valores de | forne- 
cidos pelas leis 3 e 4 (figuras 2 e 3) as quais, deverão 


2 pRaR 
B=10 


B=12 


E. DOZE 


Figura 2 — Função o =1. th). para K= 70, no caso 4. 


2. B04B 


E 2840 


E. nOZa 


[é] l Fi E | 4 E 


he th) 
Figura 3 — Função e = É, hd. para e = 2.6 mm, nos casos 


Se 6 
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Cesos 3Je5 


B-— co 


he=he (Qi) he=hu 


hezhe (0) he=he [Os] he 
Figura 4 — Andamento das curvas | = 1 (he) nos casos 3, 

[e 

4, 5 e 6. 


ser tomadas em consideração quando se adopta um 
coeficiente de segurança. 


-— À aplicação da lei 1 não é de aconselhar em 
virtude da facilidada de utilização das leis 3 e 4, sob;e- 
tudo com as actuais mini-calculadoras electrónicas. 


2,2 — Secção não rectangular 


Para canais com secção não rectangular a expres- 
são do numero de Froude é semelhante às indicadas 
em (1) e (2) com a diferença da se utilizar hm em vez 
de h, . Utilizando as leis do Quadro 2 é possível fazer- 
-S8 O traçado das funções ” o prever a evolução do 
número de Froude., Utilizando a lei 3 pode veriiicar-se 
[3] [4] o seguinte: 


Secção triangular : comportamento análogo ao caso 


3, sem valo; em 


Secção circular : comportamento análogo ao caso 
4, com valor lem para uma 


razão altura de escoamento/raio 
igual a 0.59, 


Secção parabólica : comportomento análogo ao caso 
4, com valor Im para uma de- 
terminada altura de escoamento. 


Secção trapezoidal : compo-tamento dependente da 
inclinação dos taludes e da la;- 
gura de rasto, Para um angulo 
com a normal inferior a 25º e 
uma razão entre a altura de es- 
coamento e a largura de rasto 


ITO 


inferior à unidade, o comporta- 
mento é análogo ao caso 4, 
Para um valo: da referida razão 
supe:ior à unidade, não há valor 
PeM mas sim, um máximo. Para 
angulos superiores a 25º o com- 
portamento é análogo ao caso 3. 


2.3 — Limite superior do número de Froude em 
escoamento uniforme. 


Quando existe um valor à. mínimo (im) e quan- 
do O varia entre Oie Os (figura 4), o número de Frou- 
de varia entre a unidade e um valor limite F, o qual, 
é função de i, de B, do expoente y de i na lei de resis- 
tência e do valor im Efectivamente, o número de Frou- 
de admite as seguintes transformações, válidas para 
secções iectangulares ou não rectangulares : 


U U ; 
U À 


VS Mm e (5) 


Os valores das velocidades e declives críticos U'c 
q i'c dizem respeito a um caudal fictício, ao qual corres- 
pondo uma altura crítica média de valor idêntico ao da 
altu:a uniformo média correspondente ao caudal que se 
escoa no canal, O valo: Fr co responderá, para um 
determinado i, ao caudal cuja altura h,. = om foF 
neça o mínimo i.m da função i - Nesta conformidade, 
a exp.essão do limite superior do número de Froude, 
em escoamento uniforme, é a seguinte: 


|| a 
F 
ho E ) 
cM (6) 


O valo: F, verifica-se quando se escoa um caudal 
O para o qual a respectiva altura uniforme média cumpre 
a seguinto igualdade : 


ço CM E ko (Mm) (7) 


proveniento da função E = le (Nem) correspondente à 
lei de resistência utilizada, No caso do escoamento ser 
lamina, y = 1,0 e para um escoamento puramente tur- 
Eulento, = C.5, 

Ve.ifica-se, pela presente exposição, que a expres- 
são do valor F, não depende, de uma maneira explicita, 
nem da forma da secção, nem da rugosidade de peri- 
metro molhado. Pa:a todas as formas de secção que 
edmitam a existência de valores estacionários nas fun- 
ções 1. e F, é sempre possivel, pelo menos teoricamen- 
ic, a determinaçãa fo-mal das expressões corespon- 
dentes a im ea F . O caso 4 pode considera:-se de 
fácil tratamento analítico vindo pa:a a primeira expressão 


KB ló (8) 
a qual intioduziua na expressão (2) conjuntamente com 


a lei da resistência 3, permite obter a seguinte expres- 
são de F 
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(9) 
cuja derivada, em ordem a rm , Se anula para um 


valo; desta variável igual a “ao qual substituido em 


(9) conduz à exp-essão já deduzida para F, . E de 
referir que este caso particular foi já apresentado em 
[4] não tendo, contudo, merecido da parte dos autores 
do artigo citado, a generalização que o assunto permite. 

O cálculo analítico dos valores de lem é Fj. em 
cada caso particular, não parece ter grande interesse. 
Efectivamente, sendo a expressão (6) aplicável a qual- 
que: forma da secção ou lei de resistência, a determi- 
nação de i .w deverá ser efectuada por via gráfica, em 
virtude do grande interesse prático que as curvas, do 
tipo das apresentadas, possuem e do discutivel rigor 
nas escolhas da lei a utilizar e do coeficiente de rugo- 
sidada a adoptar. A representação gráfica da função 
i permite uma rápida averiguação do regime de escoa- 
mento para um determinado O, do sentido da variação 
da F quando se altera o caudal ou as características do 
canal e o cálculo iterativo de h,., ou h., (assunto 
tratado a seguir, na secção 3). 

Em escoamento não uniforme a expressão (5) é 
válida desde que i seja substituido por J. Atendendo a 
quo a perda de carga unitária deixa de ser constante, 
deixa, igualmente, de existir um limita superior para 
o número de Froude. 


3 — DETERMINAÇÃO DA ALTURA UNIFORME 


A função i. pode ser utilizada na determinação 
da altura uniforme através das expressões que se podem 
deduzir para o número de Froude. Efectivamente, o 
parâmetro em causa, pode ter a seguinte expressão, 
para um determinado caudal O, 


E 
Q hem (6) 
É eee “À NBu 
E Ni va À 3 | 
u “um 4/8 
um Riga (10) 


a qual, em conjunto com a (5), permite a determinação, 
por via iterativa, da altura Mm é se se conhecer a função 
i-  correspondento à forma da secção. Para secções 
rectangulares, as alturas médias coincidem com as 
alturas reais e Bc é igual a Bu, tornando-se o cálculo 
relativamente fácil (em 4 apresenta-se um exemplo 
numérico). Dum modo semelhante, pode-se determinar 


Nem se se conhecer Dai 


4 — APLICAÇÕES E EXEMPLOS NUMÉRICOS 


O conhecimento do limite superior do número de 
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Froude pode ter muito interesse para o projectista, 
em virtude de muitos dimensionamentos, nos escoamen- 
tos em superfícis livre, dependerem desse parâmetro. 
Citam-se, entre outros, os seguintes problemas: condi- 
ções de entrada numa bacia de dissipação, estabilidade 
externa do escoamento e condicionalismos na variação 
de largura de um canal ao longo do seu percurso. 


O dimensionamento de um canal no qual o número 
de Froude não possa ultrapassar um determinado valor, 
torna-se possível desde que a função 1, admita um 
minimo absoluto, no intervalo de caudais de interesse. 
As figuras apresentadas, ou outras correspondentes a 
diversas formas de secção, permitem a escolha de B e 
de i mais convenientes. As mesmas figuras possibilitam 
o cálculo da altura uniforme para um determinado caudal. 


Ap:esentam-se alguns exemplos numéricos em ca- 
nais com secção rectangular com K = 70, lei de resis- 
tência 3 e escoamentos uniformes (caso 4). 


Exemplo 1 — Pretende-se que o número de Froude 
nunca exceda o valor 2. Da expres 
são (6) virá i<4 i-m Para uma 
largura B= 6 m o declive do canal 
terá que ser | < 0.011 B (figura 2). 
Para uma largura B = 4m virá para 
o declive i <. 0.0134. 


Exemplo 2 —- Num canal com i = 0.0060, O = 12 
mº si a B=4m, pretende-se deter- 
minar o valor de h,. Ao caudal indi- 
cado corresponde uma altura crítica 
he = 0.97 m. Por via iterative, ob- 
têm-se os valores = h'. = 0,80 
m e i'c = 0.00337. Pela expressão 
(5) corresponde um número de Frou- 
da de F = 1.33, Utilizando o mesmo 
valor de h, virá, pela expressão (10), 
F = 1.34. Para um cálculo expedito 
pode considerar-se como boa apro- 
aproximação h, = 0.80 m. 


Exemplo 3 — Para o mesmo canal do exemplo 2 
(com | = 0.0060 e B = 4m) preten- 
de-se determinar o valor F, . Pela 
figura 2 ou expressão (8), obtém-se 
icM = 0.00336 e, pela expressão 
(6), obtém-se F, = 1.336. 


5 — CONCLUSÕES 


Às considerações expostas permitem concluir que, 
para determinadas formas e dimensões das secções, 
existe um valor limite superior do número de Froude 
cuja expressão geral é de fácil aplicação desde que se 
conheçam as funções do declive crítico com a altura. 
Este limite superior existe quando há um valor mínimo 
para o declive crítico. À representação gráfica das refe- 
ridas funções permite uma fácil compreensão da evolu- 
ção do número de Froude e podem ser um útil auxiliar 
no dimensionamento de canais de leito fixo, em escoa- 
mento uniforme. 
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7 — NOMENCLATURA 


F  — Número de Froude; 
F,j — Limite superior do número de Froude : 


h, — Altura uniforme : 


AR — Altura uniforme média ; 


h. — Áltura crítica ; 


Nie — Altura crítica média ; 

h'.  — Altura crítica com valor igual ao da altura 
uniforme ; 

Nei — Altura crítica média com valor igual ao da 


altura uniforme média ; 
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O  —. Caudal: 

Qi — Caudal limite inferior do regime rápido; 
Os — Caudal limite superior do regime rápido : 
R — Raio hidráulico ; 

—. Declive do canal ; 

lc — Declive crítico ; 

a — Declive crítico correspondente ao caudal O; 


em — Declive crítico mínimo ; 


. — Declive crítico correspondente a he ou h' 


E cm 


B  — Largura da secção rectangular; 

B, — Largura superficial em escoamento crítico; 

Bu — Largura superficial em escoamento uni- 
forma ; 

C, — Coeficiente de Chézy, independente de R'; 

C — Coeficiente dae Chézy; 

K  — Coeficiente de Gauckler ; 

e -—— Rugosidade absoluta ; 

gq  — Aceleração da gravidade; 

U —. Velocidade média; 

U'c — Velocidade média correspondente a h'c ou 
h'em:; 

J — Perda de carga unitária ; 

et — Coeficiente de Coriolis ; 

Y  — Bcose, com 4 Coeficiente de Jaeger Man- 
zanares e 6 angulo da fundo com a hori- 
zontal. 


TÉCNICA 432 


“Do Conhecimento e da acção” 


COMENTARIO 


Respondendo aos comentários que o Eng. Medina Mar- 
tins fez o favor de fazer ao artigo «Do conhecimento e 
da acção» publicado no n.º 429 da Revista «Técnica» e 
porque esses comentários vão permitir um esclareci- 
mento mais completo do pensamento que enformou o 
referido artigo, venho solicitar, com a autorização do 
referido Engenheiro, à Direcção da «Técnicar a publi- 
cação do seguinte: 


1.º Comentário : 


A descrição dos «seres» tal como apresentados no 
artigo, afiguram-se estáticos e não dinâmicos. 


Resposta : 


Um dos «seres instrumentos» a que é necessário 
reco rer na maioria dos problemas reais é o cronômetro. 

Assim, um dos atributos do «universo dos seres» é 
o Tempo. — fcronos). 


O Cronômetro e o respectivo atributo, o tempo, 
pode ser usado de duas formas típicas : 


a) No espaço imagem do «se; em observação» 
não está incluido o tempo. 


Neste caso o tempo é o atributo de um outro «ser», 
o crometro, 

E então possivel estabelecer uma correspondência 
entre as imagens do «ser em observação» e as do crome- 
tro que desempenha aqui a função de «ser» auxiliar. 

Esta correspondência vai autorizar construir imagens 
derivadas ou integ.ais em relação ao tempo. 

O Tempo toma a configuração de um atributo uni- 
versal e independente dos objectos em observação. 

Na mecânica clássica, todos os crometros são sin- 
cronos, o tempo é universal. 


b) No espaço imagem do «ser em observação» 
está incluido o atributo — tempo. 

O «ser em observação» também é ou também inclui 
um crometro que lhe é próprio e que dá origem ao 
tempo desse ser. 

Pode até comparar-se o «tempo do sery com o 
«tempo» dado pelo cronómetro do laboratório. 
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A teoria da relatividade prevê que os dois «tempos» 
não sejam síncronos e oferece um conjunto de justifi- 
cações: acelerações, velocidades elevadas etc, 

O que importa reter na hipótese b) é a distinção 
entre o «tempo do ser» e o «tempo do cronómetro do 
laboratório». 

Resumindo a resposta ao 1.º comentário diremos 
que a discrição dos «seres» é estática ou dinamica con- 
forme o atributo tempo é tomado em linha de conta, quer 
como atributo de um ente autónome — o cronometro 
— quer como uma propriedade do próprio «ser em obser- 
vação». 


2º Comentário : 


Como intervem o conceito de «invarianciar de um 
dado «se» ? 


Resposta : 


O conceito de invariancia é sempre relativo. 

O teste da invariancia envolve as seguintes opera- 
ções : 

Definição do sistema de funções ou funcionais da 
«imagem» do «ser» cujas invariancias se deseja verificar. 

Definição de um outro sistema de funções ou fun- 
cionais da imagem do «ser» que se modifica arbitrária- 
mente, embora possam essas modificações estar sujeitas 
a um conjunto de restrições (holónomas ou anholónomas. 

Formadas estas duas classes há que testar se todos 
os membros da 1.º classe se mantêm invariantes, quendo 
um, vários ou todos os menbros da 2.º classe tomam 
valores arbritários apenas limitados pelas restrições 
holónomas ou anholónomas referidas. 

Um exemplo é o da invariancia temporal, onde o 
tempo é um membro da 2.º classe, que se faz variar 
arbitráriamente, nomeadamente invertê-lo. 

Outro exemplo, é a dos vectores-n, como o volume, 
em espaços equivoluminares; com efeito, o tensor q” é 
invariante na classe X” dos espaços equivoluminares. 

No entanto a invariancia de um atributo em rela- 
ção a um outro não significa, Obviamente, que o seja 
em relação a todos os outros; daí que a invariancia 
tenha carácter relativo. 
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Resumo dos Artigos publicados na «Técnica» n.º 432 


Ano LI — Janeiro 1976 


C. D. U. 532.516:518 


E. R. ARANTES E OLIVEIRA 
CONVERGÊNCIA DO METODO DOS ELEMENTOS FINITOS 


APLICADO AOS ESCOAMENTOS VISCOSOS 
Técnica N.º 432 — L| — 11976, p. 139 a 154 


Resumem-se primeiramente os resultados de uma teo- 
ria dos métodos  variacionais cuja utilização  per- 
mito à análisa da convergência do método dos elementos 
finitos aplicado ao escoamento de um fluído viscoso regido 
pelas equações de Navier-Stokes. O método empregado re- 
quer a decomposição do fluído em dois corpos mais sim- 
ples: um fluído viscoso sem inércia e um fluído com inércia 
mas sem viscosidade. Apresenta-se para este último um 
novo principio de minimo, 


C. D. U. 624.074 
J. 3. BRAZÃO FARINHA 
CALCULO DE ESTRUTURAS RETICULADAS NO ESPAÇO 
Técnica N.º 432 — LI — 1-1976, p. 155 a 166 


Expõe-se sucintamenta a aplicação do método de 
Cross na resolução de estruturas espaciais, constituídas 
por barras de momento de inércia constante ou variável, 
da eixo rectilineo, ortogonais entra si, ou não, 

Apresentam-sa alguns exemplos de aplicação. 
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ANTONIO BETAMIO DE ALMEIDA 


LIMITE SUPERIOR DO NUMERO DE FROUDE EM CANAIS 


PRISMÁTICOS E ESCOAMENTO UNIFORME 
Técnica N.º 4392 — L| — 1-1976, p. 167 a 172 


O autor estuda a variação do número de Froude em 


canais prismáticos e escoamento uniforme a a possibi- 
lidade de o mesmo admitir um limite superior, Apresen- 
tam-se, igualmente, algumas aplicações numéricas. 
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U. D. C. 532.54 
ANTONIO BETÂMIO DE ALMENDA 


FROUDE NUMBER UPPER BOUND IN PRISMATIC 
CHANNELS AND UNIFORM FLOW 


Técnica N.º 432 — LI —1-1976, p. 167 a 172 


The author studies Froude number range in prismatic 
channels and uniform flow and the possibility of an upper 
bound for it. Numerical aplications are also shown. 


Li-January 1976 


U. D. C. 532.516:518 


E. R. DE ARANTES E OLIVEIRA 


CONVERGENCE OF FINITE ELEMENT SOLUTIONS 
IN VISCOUS FLOW PROBLEMS 


Técnica N.º 432 -—L| -—- 14976, p. 139 a 154 


Results of a theory. variational problems are first 
summarized. Tha theory is then applied to the analisys of 
convergence of the finite element method applied to the 
flow of a viscous fluid governed by Navier-Stokes equa- 
tions. The approach requires the decomposition of the 
body into two simpler ones: a viscous non-inertial fluid 
and an inertial non-viscous one, the mechanical behaviour 
of which is governed by minimum princip'es. A new mi- 
nimum principle is presented for the inertial non-viscous 
Huid, 


U. D, C. 624,074 
J. 5. BRAZÃO FARINHA 
ANALYSIS OF SPACE FRAMES 
Técnica N.o 432 — L| — 1-1976, p. 155 a 166 


| is succintly stated the application of the Cross 
method on the resolution of rigid space structures, consti- 
tuted by bars of constant or variable inertia, of rectilinear 
axlc, orthogonal between them, or not. 

Some appíication examples are presented. 
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